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I  EXECUTIVE  SUMMARY 


Between  May  1976  and  August  1979--the  active  data-taking  phase  of 
the  Wideband  Satellite  Experiment--an  extensive  data  base  was  accumulated 
and  analyzed  to  determine  the  temporal  coherence,  spatial  coherence,  and 
frequency  coherence  of  naturally  occurring  ionospheric  scintillation  and 
its  global  morphology.  The  unique  feature  of  the  Wideband  satellite  ex¬ 
periment  was  the  S-band  phase  reference,  which  permitted  measuring  phase 
as  well  as  amplitude  scintillation  at  L  band,  UHF,  and  VHF. 

For  striations  greater  than  ~200  m,  where  diffraction  effects  can  be 
safely  neglected,  the  Wideband  phase-scintillation  structure  is  an  inte¬ 
gral  measure  of  the  in-situ  structure  along  the  propagation  path.  Ii 
the  shape  of  the  in-situ  spectrum  of  the  irregularities  does  not  vary 
along  the  propagation  path,  then  the  power-law  index  of  the  phase  spectrum 

is  one  plus  the  corresponding  index  of  the  one-dimensional  in-situ  spec- 

-2  -3 

tral  index.  Thus,  a  k  in-situ  spectrum  corresponds  to  a  k  phase 

spectrum. 

The  Wideband  satellite  data,  however,  have  consistently  shown  shal¬ 
lower  spectral  indices  than  expected  that  vary  systematically  with  chang¬ 
ing  perturbation  strength.  Only  recently,  in-situ  measurements  from  the 
Atmospheric  Explorer-E  (AE-E)  satellite  and  the  PLUMEX  experiment  have 
demonstrated  that  the  Wideband  data  correctly  predicted  the  in-situ 
structure.  The  PLUMEX  rocket  data,  moreover,  have  extended  the  measure¬ 
ments  to  structures  well  below  200  m.  There,  the  in-situ  spectra  roll 
off  more  rapidly  than  expected. 

The  in-situ  spectral-density  function  evidently  has  a  two-component, 
power-law  shape.  For  structures  smaller  than  ~200  m,  where  striation 
shape  dominates,  the  spectra  fall  off  more  rapidly  than  k  .  For  struc¬ 
tures  larger  than  ~200  m,  where  the  striation-size  distribution  dominates, 

"2 

the  spectra  fall  off  less  rapidly  than  k  .  This  characterization  strictly 
applies  only  to  the  most  intensely  disturbed  (high-electron-density) 
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regions.  Bottomside  spread  F  and  density  fluctuations  within  equatorial 

-2.3  ±  0.2 

depletions  show  the  expected  k  *  '  power-law  behavior,  but  do  not 

contribute  significantly  to  disturbances  on  propagation  paths  that  also 
intercept  highly  structured,  high-density  regions. 

The  Wideband  satellite  data  are  clearly  dominated  by  the  large-scale 
structures,  which  are  generally  larger  than  the  Fresnel  radius.  This 
effect  has  been  intuitively  clear  from  the  earliest  Wideband  measurements. 
Under  this  contract,  we  have  carefully  analyzed  the  mutual  coherence 
function,  which,  in  principle,  is  unaffected  by  the  diffraction  effects 
that  distort  the  phase  data.  The  results  show  clear  evidence  of  a  sys¬ 
tematic  steepening  of  the  integrated  phase  spectrum  at  high  frequencies 
that  is  consistent  with  the  PLUMEX  data. 

We  have  also  analyzed  the  total  electron  content  (TEC)  data  recorded 
at  Poker  Flat,  Alaska  and  the  equatorial  stations  at  Kwajalein,  Marshall 
Islands,  and  Ancon,  Peru.  The  analysis  of  the  auroral-zone  data  was  con¬ 
centrated  on  the  occurrence  and  structure  of  distinct  (TEC)  boundaries 
that  are  evidently  associated  with  an  enhanced,  unstable  F  region.  The 
scintillation  signature  of  these  regions  shows  a  unique  sheet-like  aniso¬ 
tropy  and  recently  has  generated  considerable  interest.  The  equatorial 
TEC  patterns  are  much  more  systematic,  and  reveal  a  considerable  amount 
of  information  about  the  distribution  of  equatorial  ionization  and  its 
seasonal  and  solar  cycle  variation. 

Both  the  propagation  effects  and  the  background  ionization  structure 
have  ramifications  for  predicting  propagation  effects  in  disturbed  nuclear 
environments.  The  two-component  spectral  structure  illustrates  the  im¬ 
portance  of  establishing  where  the  spectrum  becomes  distribution  dominated 
rather  than  edge  dominated.  The  former  regime  dominates  the  Wideband  scin¬ 
tillation  structure,  whereas  the  latter  regime  is  most  readily  amenable 
to  numerical  simulation  and  theoretical  analysis.  The  analysis  of  the 
background  ionization  structure  is  important  for  establishing  the  rela¬ 
tionships  between  the  macroscale  and  the  various  microscale  structure 
regimes. 
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II  INTRODUCTION 


The  active  data-taking  phase  of  the  Wideband  Satellite  Experiment 
was  formally  completed  in  August  1979  [Rino  et  al , ,  1980].  Extensive 
data  were  accumulated  and  analyzed  to  determine  the  temporal  coherence, 
spatial  coherence,  and  frequency  coherence  of  naturally  occurring  scin¬ 
tillation  and  its  global  morphology.  The  systematic  analysis  of  these 
data  has  provided  a  consistent  characterization  of  the  underlying  irregu¬ 
larity  structure  to  scale  sizes  as  small  as  a  few  hundred  meters,  beyond 
which  the  effects  of  diffraction  severely  distort  the  phase  data.  The 
PLUMEX  rocket  experiment  has,  moreover,  provided  in-situ  measurements  to 
scale  sizes  approaching  one  meter  [McDaniel,  19801.  Many  of  the  results 
were  unexpected,  however,  and  it  is  useful  to  review  the  theory  that  has 
guided  the  interpretation  of  these  experimental  data. 

It  is  now  generally  accepted  that  convective  instabilities  cause  the 
development  of  an  enhanced  continuum  of  intermediate-scale  plasma  struc¬ 
tures.  The  gradient-drift  and  Rayleigh-Tay lor  instabilities  are  well- 
known  examples  of  convective  instabilities.  When  a  convective  instability 
is  operative,  an  initially  unstable  gradient  first  steepens  and  then 
finger-like  striations  or  depletions  develop.  The  size  of  the  primary 
striations  or  depletions  can  vary  considerably. 

An  important  finding  from  numerical  simulations  of  barium  striation 
development  is  that  the  size  of  the  primary  striations  is  controlled  by 
the  initial  gradient  [Keskinen  et  al , ,  1979],  When  the  unstable  plasma 
is  completely  isolated,  the  primary  striations  bifurcate,  thereby  generat¬ 
ing  smaller  structures  which  themselves  bifurcate.  This  cascade  process 
continues  until  the  fine  structure  is  small  enough  to  be  rapidly  eroded 
by  classical  cross-field  diffusion. 

The  main  spectral  characteristics  of  the  structure  continuum  can 
readily  be  inferred  when  such  a  process  is  operative.  Photographic  data 
suggest  that  the  striation  cross-section  is  essentially  rectangular  ("top-hat" 
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_2 

like),  which  leads  directly  to  a  k  one-dimensional  spectral-density  func¬ 
tion  (SDF)  wherever  the  steep  edges  dominate.  The  canonical  SDF  is  then 

_2 

k  between  an  outer-scale  wave  number  corresponding  to  the  primary  stria- 
tion  size  and  an  inner-scale  wave  number  corresponding  to  the  ion  gyrora- 
dius. 

An  isolated  barium  cloud  does  not  conform  rigorously  to  this  simple 

picture,  however.  Numerical  simulations  by  Keskinen  et  al.  [1981]  have 

_2 

shown,  for  example,  that  a  k  SDF  actually  occurs  only  for  SDFs  measured 
across  the  tips  of  the  primary  striation.  The  cross  gradient  SDFs  are 
somewhat  steeper  with  a  spectral  index  between  2.2  and  2.5.  The  steeper 
slopes  persist  when  the  structure  achieves  a  quasi-steady  isotropic  state. 
These  spectral  characteristics  agree  well  with  bottomside  spread  F  mea¬ 
sured  during  the  PLUMEX  experiment  [Keskinen  et  al.,  1981], 

Rapid  progress  is  also  being  made  in  understanding  the  various  mecha¬ 
nisms  that  generate  the  meter-to-submeter  structure  responsible  for  en¬ 
hanced  radar  field-aligned  backscatter  echoes  [Ossakow,  1981],  The  spec¬ 
tral  continuum  for  naturally  spread  F,  however,  is  more  complex  than  the 
spectral  continuum  for  an  isolated  barium  cloud.  For  example,  the  power- 
law  spectral  indices  deduced  from  the  Wideband  satellite  data  have  been 
consistently  less  than  2,  not  greater  than  2  as  would  be  expected  from 
the  totality  of  numerical  simulation  studies. 

The  estimates  of  the  in-situ  spectral  index  are  derived  from  the 
phase-scintillation  spectra  and  are  consistent  with  the  behavior  of  the 
intensity  scintillation  and  frequency  coherence  of  the  signal  under 
strong  scatter  conditions.  Nonetheless,  the  scintillation  structure  re¬ 
sponds  to  the  integral  of  the  structure  along  the  entire  propagation 
path.  Because  of  this  "geometrical  smearing,"  it  was  initially  thought 
that  the  scintillation  structure  gave  a  biased  estimate  of  the  in-situ 
spectral  index. 

Recent  comparisons  of  Wideband  satellite  data  with  nearly  coincident 
in-situ  measurements  using  the  AE-E  satellite  have  shown,  however,  that 
the  indirectly  inferred  shallow  and  systematically  varying  spectral  in¬ 
dices  are  indeed  characteristic  of  the  in-situ  irregularity  structures 
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with  scale  sizes  between  several  hundred  kilometers  and  several  hundred 
meters  [Livingston  et  al.,  1981],  The  measured  and  inferred  perturbation 
strengths  show  systematic  differences,  which  are  attributed  to  the  fact 
that  the  AE-E  satellite  sampled  only  a  single  altitude.  This  feature  was 
demonstrated  by  the  phase  and  in-situ  data  from  the  PLUMEX  I  rocket,  which 
showed  nearly  perfect  agreement  for  a  near  vertical  path  [Rino  et  al. , 
1981]. 

To  reconcile  the  intermediate-scale  measurements,  which  show  a  much 
more  shallowly  sloped  SDF,  with  the  more  steeply  sloped  transition-scale 
structures,  we  must  postulate  a  discontinuity  in  the  power-law  index  be¬ 
tween  the  spatial  frequencies  corresponding  to  ~500  m  and  ~200  m.  The 
PLUMEX  I  data  have  provided  direct  evidence  of  such  a  two-component  spec¬ 
trum  [Rino  et  al. ,  1981;  Kelley  et  al. ,  1981],  Even  the  transitio-scale 
structures,  however,  show  a  variable  spectral  index,  in  some  regions 
steepening  to  values  as  high  as  4. 

Such  departures  from  the  expected  behavior  of  isolated  convective 
instabilities  with  uniform  driving  conditions  are  evidently  caused  by 
the  very  complex  environment  in  which  the  local  instability  operates,  as 
has  been  suggested  by  Wittwer  [l981_j.  With  regard  to  the  transition- 
scale  structure,  evidence  mounts  that  the  enhanced  diffusion  inhibits  the 
bifurcation  process  so  that  the  primary  structures  persist  much  longer 
than  expected--the  "freezing"  phenomenon  [McDonald  et  al..  19 80 J . 

The  cause  of  the  anomalous  diffusion  is  believed  to  be  E-region 
coupling,  which  undoubtedly  plays  an  important  role  in  establishing  the 
structure  of  auroral-zone  scintillation  _ Vickrey ,  1981_.  The  hypothesis 
seems  reasonable  that  anomalously  enhanced  diffusion  acting  to  scale 
sizes  approaching  100  m  or  more  can  explain  the  observed  spectral  steep¬ 
ening  over  that  measured  in  isolated  bottom,  side  structures. 

Turning  to  the  intermediate-scale  structures,  we  first  note  that  in 
the  PLUMEX  I  data  they  were  most  intense  in  the  high-density  regions  ad¬ 
jacent  to  the  large-scale  depletions  that  are  delineated  by  enhanced 
radar  backscatter.  Indeed,  evidence  mounts  that  the  west  wall  of  the 
depletion  is  convcctively  unstable  because  of  an  eastward  neutral  wind 
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[Tsunoda,  1981;  Singh  et  al. ,  1981],  Because  the  gradient  varies  with 
height,  the  locally  unstable  regions  should  have  a  varying  local  outer 
scale.  Wittwer  [ 1981 J  argues  that  the  superposition  of  the  spectra  from 
many  such  regions  can  account  for  the  shallow  and  varying  spectral  index 
that  characterizes  the  intermediate-scale  irregularity  structures. 

To  summarize,  therefore,  the  Wideband  satellite  data,  the  more  recent 
PLUMEX  data,  and  the  structure  of  isolated  convective  instabilities  can 
be  reconciled,  at  least  qualitatively.  The  break  in  the  spectral  slope 
is  formally  the  point  beyond  which  the  SDF  becomes  dominated  by  the  stria- 
tion  profile  and  the  effects  of  anomalously  enhanced  diffusion  manifests 
itself.  In  the  shallowly  sloped  region,  the  SDF  shape  is  determined  by 
the  configuration  of  the  locally  striated  subregions.  As  yet,  however, 
these  characteristics  have  not  been  verified  by  simulations. 

The  potential  ramifications  for  systems  are  important  because  the 
scintillation  structure  under  strong-scatter  conditions  is  a  sensitive 
function  of  the  power-law  index.  Scintillation  effects  are  generally 
caused  by  structures  comparable  to  the  Fresnel  radius  and  smaller.  This 
is  strictly  true,  however,  only  for  intensity  scintillation  under  weak- 
scatter  conditions.  Under  strong-scatter  conditions,  structures  larger 
than  a  Fresnel  radius,  but  smaller  than  the  outer  scale,  dominate  the 
scintillation  effects;  moreover,  the  spectral  characteristics  of  the 
underlying  irregularity  structure  manifests  itself  most  clearly  under 
strong-scatter  conditions. 

The  Wideband  satellite  data  behave  as  if  only  the  large-scale, 
shallowly  sloped  spectral  regime  were  present.  This  behavior,  evidently, 
is  due  to  the  transition  region  scale  sizes  in  the  transition  region  are 
larger  than  the  Fresnel  radius  only  at  L  band  and  above,  where  saturated 
scintillation  conditions  are  rarely  observed.  Under  this  contract,  we 
have  pursued  a  more  refined  analysis  of  the  disturbed  Wideband  satellite 
passes  to  verify  these  concepts. 

In  a  paper  submitted  for  publication  by  C.  L.  Rino  and  C.  H.  Liu 
(1981),  we  develop  a  simple  set  of  parameters  for  the  scintillation  index 
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that  can  be  used  to  determine  the  characteristics  of  the  dominant  spectral 
regime  when  only  intensity-scintillation  data  are  available.  The  results 
reconfirm  the  shallow  slopes.  Midlatitude  data  from  ATS-6  recorder  at 
140  MHz,  however,  show  different  characteristics.  This  is  most  likely 
caused  by  an  outer-scale  cutoff,  but  possibly  caused  by  spectral  steepen¬ 
ing.  In  any  case,  it  underscores  the  variability  of  naturally  occurring 
scintillation  producing  irregularities. 

To  test  the  Wideband  satellite  data  for  some  manifestation  of  the 
transition-scale  regime,  we  have  also  analyzed  the  mutual  coherence  func¬ 
tion.  The  mutual  coherence  function  is  important  in  its  own  right,  be¬ 
cause  in  the  Rayleigh  scatter  limit,  it  completely  characterizes  the 
fading  structure  at  a  single  frequency.  In  theory,  moreover,  the  mutual 
coherence  function  is  independent  of  the  Fresnel  radius.  Thus,  its  func¬ 
tional  form  should  directly  relate  to  the  underlying  irregularity  struc¬ 
ture. 

Our  results  show  that  this  is  strictly  true  only  under  conditions 
of  strong  scatter.  Under  weak-scatter  conditions  the  slow,  trend-like 
phase  variations  distort  the  mutual  coherence  function.  The  results  have, 
nonetheless,  provided  evidence  of  a  break  in  the  spectral  slope  that  can¬ 
not  be  detected  in  the  phase  data  alone  because  of  diffraction  effects. 

In  Section  III,  the  results  of  an  analysis  of  the  total  electron  con¬ 
tent  (TEC)  data  from  the  Poker  Flat  station  are  reported.  We  have  concen¬ 
trated  on  the  occurrence  and  structure  of  distinct  TEC-boundary  features 
to  determine  their  association  with  the  localized  high-latitude  scintilla¬ 
tion  feature  that  is  currently  of  considerable  interest  [Rino,  1980; 
Vickrey,  1981]. 

Distinct  TEC-boundary  features  occur  in  approximately  50  percent  of 
the  nighttime  passes,  and  are  generally  manifestations  of  an  enhanced,  un¬ 
stable  F  region.  The  TEC  gradient  is  positively  correlated  with  magnetic 
activity  during  the  summer  months,  but  the  correlation  reverses  sense 
during  the  winter  months.  The  occurrence  of  phase-scintillation  enhance¬ 
ments  within  the  boundary  region  is  positively  correlated  with  the  magni¬ 
tude  of  the  TEC  gradient. 
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Section  IV  reviews  the  results  of  TEC  studies  at  two  equatorial  sta¬ 
tions,  Ancon,  Peru  in  the  American  sector  and  Kwajalein,  Marshall  Islands 
in  the  Pacific  sector.  A  consistent  F-layer  model  is  established  from 
which  the  long-term  solar  cycle  and  seasonal  dependences  of  TEC  can  be 
analyzed.  The  overall,  large-scale  electrodynamic  behavior  at  both  sta¬ 
tions  is  consistent  with  an  F-region  dynamo  driven  by  thermospheric  winds 
The  dynamo  action  enhances  vertical  ionization  drift  in  the  evening  and 
generates  the  irregular  latitudinal  distribution  of  TEC  known  as  the  equa 
torial  anomaly.  During  the  equinoxes,  the  TEC  shows  a  strong  solar-cycle 
dependence.  Also  during  these  periods,  the  effects  of  meridional  neutral 
winds  on  the  anomaly  symmetry  are  clear.  During  the  solstices,  when  the 
dynamo  is  weak,  the  solar  effects  are  less  dominant,  and  local  influences 
seem  to  control  TEC  behavior.  The  association  between  TEC  and  scintilla¬ 
tion  is  also  considered.  As  expected,  enhanced  TEC  is  necessary  for,  but 
does  not  ensure,  scintillation  occurrence.  Detailed  comparisons  suggest 
that  the  scintillation  occurs  when  the  large-scale  TEC  distributions  are 
changing. 
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III  THE  MUTUAL  COHERENCE  FUNCTION 


A.  Introduction 

Under  conditions  of  narrow-angle  scattering,  a  wave  field,  U(r), 
that  interacts  with  a  randomly  irregular  medium  admits  the  representation 

U(r)  =  Re[u(7,  z)  exp{ikz}]  ,  (1) 

where  u(p,  z)  satisfies  the  parabolic  wave  equation.  In  Eq.  (1),  the  z 
axis  lies  along  the  propagation  direction.  The  mutual  coherence  function 
is  generally  defined  as 

Ma7;  z)  4  <u(7,  z)  u*  (£'  z)>  ,  (2) 

because  the  correlation  of  the  wavefield  along  z  can  be  neglected  for 
most  applications. 

We  have  tacitly  assumed  that  is  independent  of  location  in  the 
transverse  plane,  although  this  is  never  strictly  true  in  practice.  In¬ 
deed,  under  the  narrow-angle  scatter  approximation,  R^  must  satisfy  the 
differential  equation 


dR 

_ u 

dz 


-i  k/2 


R  -  1/2 
u 


dD(p,  p') 
dz 


In  Eq.  (3),  D(p,  p  )  is  the  phase-structure  function 


(3) 


0(7,  7')  =  <[50(7)  -  60(7,)}2> 


(4) 


where 


6<0(P) 


N  (?,  11)  dT\ 
e 


(5) 


r  is  Che  classical  electron  radius,  X  =  2it/k  is  the  wavelength,  and 
e 

Ne(p,  z)  is  the  local  electron  density. 

If  D  is  a  function  only  of  AP»  then  the  transverse  Laplacian  opera¬ 
tors  in  Eq.  (3)  cancel  and 

Ru(&?;  z)  =  exp{-  1/2  D(AP) }  .  (6) 

The  simple  form  of  Eq.  (6)  was  known  long  before  a  full  multiple-scatter 
theory  was  developed.  Booker,  Ratcliffe,  and  Shinn  [1950]  pointed  out 
that  the  mutual  coherence  function  (MCF)  is  invariant  under  conditions 
of  free-space  propagation.  Fejer  [1953]  seems  to  have  first  derived 
Eq.  (6).  Bramley  [1954]  later  pointed  out  that  Eq.  (6)  can  be  deduced 
immediately  from  the  Gaussian  phase-screen  model. 

The  two-dimensional  Fourier  transform  of  Eq.  (6)  gives  the  spatial 
spectral-density  function  of  the  wave  field--the  so-called  angular  spec¬ 
trum.  Under  the  usual  assumption  of  "frozen"  irregularity  structures, 
moreover,  the  MCF  for  the  complex  signal  can  be  obtained  by  replacing 
AP  by  vAt,  where  v  is  the  irregularity-drift  velocity.  The  one-dimensional 
Fourier  transform  of  the  temporal  MCF  gives  the  Doppler  spectrum. 

The  most  general  form  of  the  phase-structure  function  is  too  cumbrous 
for  use  in  most  computations.  Rino  [1980],  however,  shows  that 

D(y)  ~  c6*|yfitl(2v”1’  2)  ,  (7) 

where  v  is  the  power-law  index  parameter  for  a  three-dimensional  spectral- 
density  function  of  the  form  c  q"^v+^,  The  phase-structure  constant 
can  be  derived  from  the  formula 
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(8) 


c  2  =  — £■ 

6m  2k 


(v  -  4.5v  +  5.5)  for  1  s  v  ^  2 


where 

2  2 

C  =  r  n  C 
p  e  p  s 


(9) 


In  Eq.  (9),  l  is  the  length  of  the  propa.&au 


X - - 


The  effects  of  anisotropy  can  be  accommodated  by  replacing  y  for 
v  and  multiplying  C  by  G,  where  ve££  is  the  effective  velocity 

parameter,  and  G  is  a  coherent  amplification  factor  defined  in  Rino  (1979a). 
The  main  point  is  that  a  detailed  measurement  of  the  MCF  can  provide  in¬ 
formation  about  the  underlying  irregularity  structure  that  is  unencumbered 
by  a  complicated  dependence  on  the  Fresnel  radius,  as  is  the  case  when 
intensity  or  phase  statistics  are  used. 


The  Wideband  satellite  data  provide  a  unique  opportunity  to  measure 
the  mutual  coherence  function.  The  S-band  reference  frequency  is  essen¬ 
tially  undisturbed  even  under  conditions  of  very  strong  scatter  at  VHF 
and  UHF.  Moreover,  under  saturated-scatter  conditions  (S^  ~  1),  it  is 
generally  assumed  that  Rayleigh  scattering  prevails.  In  that  case,  the 
MCF  is  sufficient  to  completely  specify  the  signal  statistics  at  a  single 
frequency. 


B.  Data  Analysis 

To  compute  the  MCF,  we  have  used  Wideband  satellite  data  sampled  at 
a  250-Hz  rate.  The  amplitude  and  phase  data  were  separately  detrended 
using  a  25-s  (0.04-Hz)  cutoff.  This  step  is  particularly  important  in 
computing  the  MCF  because  it  effectively  establishes  the  local-phase 
reference.  To  compute  the  actual  MCF,  we  first  compute  the  Doppler  spec¬ 
trum  using  a  standard  FFT  algorithm.  A  second  Fourier  transform  then 
gives  the  mutual  coherence  function  estimate. 

The  Doppler  spectrum  is  very  nearly  symmetric.  Thus,  to  force  a 
purely  real  MCF,  we  have  averaged  the  positive  and  negative  frequencies 
together  before  computing  the  final  FFT.  A  standard  window  function  was 
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also  applied  to  the  detrended  data  to  minimize  sidelobe  effects.  Examples 
of  the  Doppler  spectra  and  the  corresponding  MCFs  for  typical  moderate 
and  strong  disturbances  are  shown  in  Figures  1(a)  and  (b),  respectively. 

The  left-most  frames  show  the  Doppler  spectra  on  log-linear  plots. 

The  data  have  been  normalized  to  unit  average  intensity  over  the  16-s 
data  interval.  Thus,  the  integral  of  the  Doppler  spectrum  is  unity.  The 
reduction  in  intensity  of  the  low-frequency  spectral  components  and  the 
attendant  spectral  broadening  induced  by  strong  scattering  is  clearly 
evident.  The  middle  frame  shows  the  same  Doppler  spectra  on  a  log-log 
plo<-  For  clarity  of  presentation,  however,  the  negative  frequencies 
have  been  offset  by  20  dB. 

As  might  be  expected,  the  higher  frequencies  show  a  power-law  regime 
with  the  same  power-law  index  as  the  integrated-phase  spectrum.  The 
lower  frequencies  are  distorted,  however,  and  bear  no  simple  relation  to 
the  integrated-phase  spectrum.  Nonetheless,  the  shift  of  the  spectral 
content  to  higher  frequencies  with  increasing  perturbation  strength  is 
conspicuous  in  both  the  semi-log  and  log-log  displays. 

Now,  to  the  extent  that  Eq.  (6)  is  strictly  applicable,  the  phase 
spectrum  can  be  extracted  by  first  taking  the  logarithm  of  the  MCF  and 
then  computing  the  Fourier  transform  of  the  result.  The  problem  is  that 
the  measured  MCFs  are  not  strictly  positive  as  the  righthand  frames  of 
Figures  1(a)  and  1(b)  show.  Indeed,  for  the  weak  perturbation  the  nega¬ 
tive  excursion  appears  to  be  systematic. 

In  retrospect,  this  result  might  have  been  expected.  The  phase- 
structure  function,  in  fact,  any  measured-structure  function,  is  a  func¬ 
tion  of  AP  only  to  some  maximum  separation.  In  a  simple  homogeneous  en¬ 
vironment,  the  phase-structure  function  converges  to  twice  the  rms  phase 
at  this  point.  In  practice,  trend-like  phase  variations  upset  this  con¬ 
vergence,  and  the  structure  function  does  not  remain  position-  or  time- 
invariant  for  large  separations.  Under  such  conditions,  Eq.  (3)  implies 
a  weak  dependence  on  the  Fresnel  radius.  Under  strong-scatter  conditions, 
decorrelation  occurs  sufficiently  rapidly  that  the  non-stationary  compo¬ 
nent  can  evidently  be  neglected. 
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Because  a  single  MCF  estimate  is  itself  a  random  function,  we  must 
average  many  estimates  together  to  determine  the  functional  form  of  the 
underlying  MCF.  Each  measured  MCF  is  first  corrected  for  wavefront  curva¬ 
ture  effects  as  described  in  Rino  and  Owen  [1980(b)].  All  other  geometri¬ 
cal  effects  scale  in  exactly  the  same  way  as  does  the  phase  turbulent 
strength  parameter,  T,  described  in  Fremouw  et  al.  [1978],  [see  also  Rino 
and  Owen  (1980a)].  Thus,  the  corrected  MCFs  with  common  T  values  should 
be  statistically  similar. 

We  have,  therefore,  sorted  the  corrected  MCF  estimates  by  turbulent 
strength  as  deduced  from  the  corresponding  phase  spectra  and  averaged  the 
results.  The  upper  frames  in  Figure  2  show  the  results  derived  from  the 
more  disturbed  portions  of  24  passes  recorded  at  Ancon,  Peru.  To  the  ex¬ 
tent  that  the  form  implied  by  Eq.  (6)  and  (7)  is  correct,  a  log-log  plot 
of  the  base  e  logarithm  of  the  MCF  should  be  linear.  Thus,  in  the  lower 
frames  in  Figure  2,  we  have  generated  such  a  plot  for  the  strictly  posi¬ 
tive  portion  of  the  measured  and  averaged  MCFs,  together  with  a  linear 
least-squares  fit  to  the  resultant  curve. 

The  corresponding  weak-scatter  data  are  shown  in  Figure  3.  Here  we 
see  a  distinct  and  systematic  negative  tail.  Nonetheless,  the  strictly 
positive  portion  adheres  quite  well  to  the  exponential  power-law  model. 

The  varying  termination  of  the  MCFs  is  caused  by  the  spherical  wavefront 
correction. 

Data  recorded  at  the  Wideband  station  in  Kwajalein  have  also  been 
processed  in  this  manner.  To  summarize  the  results,  in  Figures  4  and  5, 
we  have  plotted  the  index  parameter  derived  from  the  linear  fit  to  the 
log-log  plots  of  the  base  e  logarithm  of  the  measured  MCFs.  The  time  to 
50-percent  decorrelation,  t  ,  is  also  shown.  The  straight  line  superim¬ 
posed  on  the  plot  of  v  versus  T  is  the  empirically  derived  relation  ob¬ 
tained  from  the  Wideband  AE-E  data  comparisons.  The  solid  curve  superim¬ 
posed  on  the  versus  T  plot  is  the  theoretical  curve  derived  from  Eq.  (6) 
and  Eq.  (7),  but  using  the  v  value  computed  from  the  empirical  v  versus  T 
relation  rather  than  the  measured  value. 
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2  STRONG  SCATTER  MCFs  FROM  ANCON  DATA  TOGETHER  WITH  LOG-LOG 
PLOTS  OF  ESTIMATES  OF  PHASE-STRUCTURE  FUNCTION.  The  straight  lines 
are  linear  least-squares  fits  to  the  data. 
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5  SUMMARY  PLOT  FOR  ANCON  DATA  SHOWING  TIME  LAG  TO  50-PERCENT 
DECORRELATION  AND  SPECTRAL-INDEX  PARAMETER  DERIVED  FROM  LOG 
LINEAR  LEAST  SQUARES  FIT  TO  ESTIMATE  OF  PHASE-STRUCTURE  FUNCTION 
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The  empirically  derived  v  versus  T  curve  fits  the  weakLy  to  moderately 
disturbed  Kwajalein  data  [T  <  -15  dB  in  Figure  4]  very  well.  Thus,  in 
this  range,  the  MCF  is  consistent  with  the  independently  processed  phase 
data  and  the  AE-E  satellite  in-situ  data.  For  the  more  disturbed  data, 
the  inferred  spectral  index  steepens.  Intuitively,  the  behavior  is  ex¬ 
pected  if  the  in-situ  SDF  had  an  abrupt  change  in  the  power-law  slope. 

This  feature  evidently  does  not  manifest  itself  in  the  behavior  of 
T  »  since  the  versus  T  data  follow  the  theoretical  computations  based 
on  the  empirically  derived  v  versus  T  curve,  except  for  the  last  data 
point.  The  resolution  limit  for  log^T^  is  -2.4.  In  any  case,  t  is  de¬ 
rived  from  a  single  point,  whereas  v  is  derived  from  a  range  of  MCF  values. 
Moreover,  the  behavior  of  the  MCF  for  a  two-component  SDF  need  not  give 
rise  to  a  simple  disjointed  two-component  MCF. 

In  the  Ancon  data  shown  in  Figure  5,  the  behavior  of  versus  T  is 
virtually  identical  to  the  Kwajalein  data.  However,  the  versus  T 
curve  shows  pronounced  departures  from  the  empirically  derived  T  versus 
v  curve.  We  conclude  that  when  comparatively  simple  measures  are  used, 
e.g.,  t  ,  that  MCF  theory  and  data  agree  well.  A  more  refined  measure  of 
the  MCF  structure,  however,  shows  evidence  of  a  more  complex  high-frequency 
SDF  behavior,  in  particular,  a  steepening  of  the  spectral  index. 

To  illustrate  this  effect  directly  in  Figure  6,  we  show  a  single 
power-law  SDF  superimposed  on  a  two-component  power-law  SDF  together  with 
the  corresponding  MCFs  and  the  log-log  display  of  the  corresponding  phase 
structure  function.  The  computations  were  all  done  numerically.  The 
measured  parameter  agrees  to  two  decimal  places  for  the  two  SDFs.  In 
the  log-log  display,  however,  a  steepening  of  the  power-law  slope  at 
short  lags  is  clearly  evident.  Such  behavior  completely  explains  the 
data  shown  in  Figures  4  and  5. 

As  noted  earlier,  a  Rayleigh  process  is  completely  characterized  by 
the  MCF.  In  particular,  =  1  and 

<Il'>  -  1  =  |Ru(At) | 2  ,  (10) 
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FIGURE  6  NUMERICAL  COMPUTATION  SHOWING  THE  EFFECT  OF  A  DISCONTINUITY  IN  THE  POWER-LAW  INDEX 


where  I  is  Che  signaL  intensity  normalized  to  unity  average.  Strictly 
speaking,  Eq.  (10)  is  a  necessary  but  not  sufficient  condition.  Rino 
(1979b)  shows  that  Eq.  (10)  is  the  strong-scatter  limit  as  long  as  the 
power-law  index  for  phase  is  less  than  3. 

Because  the  dominant  irregularity  structures  fall  well  within  this 
limitation,  Eq.  (10)  should  apply  to  saturated  ionospheric  scintillation. 
This  was  verified  indirectly  in  Rino  and  Owen  (1980a,  1981)  by  plotting 
t  versus  T,  where  t  is  the  time  interval  to  50-percent  decorrelation 
of  intensity.  The  data  agreed  very  well  with  theoretical  predictions 
based  on  Eqs.  (6),  (7),  and  (10). 

The  model  implies  a  simple  relationship  between  and  under 

strong-scatter  conditions.  To  verify  this,  we  have  averaged  the  corrected 

T  measurements  in  exactly  the  same  way  as  the  T  measurements  that  were 
v  I 

reported  in  Rino  and  Owen  (1981a)  were  processed.  The  results  are  shown 
in  Figures  7  and  8.  As  before,  the  solid  curves  are  derived  from  theory, 
but  use  the  empirically  derived  relation  between  v  and  T.  As  predicted, 
the  data  are  simple  displacements  of  the  corresponding  curves  in  Rino 
and  Owen  (1981). 

C.  Discussion 

In  this  section  we  have  presented  the  first  detailed  MCF  computations 
using  ionospheric  radiowave  scintillation  data.  It  has  been  known  from 
the  time  of  the  early  development  of  radiowave  scintillation  theory  that 
the  MCF  is  invariant  under  free  space  propagation  and  independent  of  the 
Fresnel  radius.  Thus,  it  has  attractive  properties  for  diagnostic  pur¬ 
poses.  Our  results  show,  however,  that  slow,  trend-like  phase  variations 
disrupt  the  simple  form  of  the  MCF,  particularly  under  weak-scatter  con¬ 
ditions.  For  shorter  time  delays  where  the  measured  MCFs  are  strictly 
positive,  however,  the  theory  fits  very  well. 

The  structure  of  the  MCF  for  short  time  lags  verifies  the  dominant 
spectral  characteristics  of  the  intermediate-scale  irregularities  inferred 
independently  from  the  Wideband  phase  data.  Moreover,  the  MCF  structure 
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FIGURE  7  PLOT  OF  TIME  LAG  TO  50-PERCENT  DECORRELATION  vs  SPECTRAL 
STRENGTH  FROM  KWAJALE1N  DATA 

showed  evidence  of  spectral  steepening  at  high  spatial  frequencies.  This 
feature  has  not  been  detected  in  the  Wideband  phase  data,  evidently  be¬ 
cause  of  diffraction  effects,  A  steepening  of  the  power-law  slope  near 
500  m  was,  however,  directly  measured  in  the  PLUMEX  I  rocket  data. 

Finally,  our  results  verify  the  use  of  the  Rayleigh  model  for  predict 
ing  the  coherence  time  of  radiowave  scintillation  under  saturated  (S^  ~  1) 
conditions.  A  straightforward  extension  of  the  model  can  be  used  to  com¬ 
pute  the  Doppler  spectrum  or  the  angular  spectrum  of  plane  waves. 
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FIGURE  8  PLOT  OF  TIME  LAG  TO  50-PERCENT  DECORRELATION  vs  SPECTRAL 
STRENGTH  FROM  ANCON  DATA 
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IV  AURORAL- ZONE  TEC  STRUCTURE 


A.  Introduction 

In  this  section  we  shall  report  the  results  of  a  summary  analysis 
of  the  high-latitude  total  electron  content  (TEC)  structure  as  measured 
by  the  Wideband  satellite.  The  scintillation  morphology  as  measured  by 
Wideband  has  received  considerable  attention,  particularly  the  prominant 
nighttime  enhancement  that  has  been  attributed  to  an  unstable  F  region 
containing  sheet-like  irregularities  [see  Rino  et  al.  (1980)  and  the 
references  cited  therein].  We  have,  therefore,  concentrated  on  the  TEC 
structures  that  can  be  associated  with  this  phenomenon.  Indeed,  the 
F-region  source  was  first  identified  by  using  multiple-station  TEC  mea¬ 
surements  [Rino  and  Owen,  1980b]. 

One  component  of  the  nighttime  auroral-zone  particle  precipitation 
sustains  an  enhanced  E  layer  over  a  broad  latitudinal  region.  The  opti¬ 
cal  signature  of  this  precipitation  is  called  the  diffuse  aurora.  The 
more  dynamic  discrete  auroras  generally  lie  poleward  of  the  diffuse  auroras, 
but  they  can  penetrate  into  the  region  of  the  diffuse  auroras  during  sub¬ 
storms  . 

The  equatorial  edge  of  the  diffuse  aurora  is  a  well-defined  feature 
in  both  optical  and  particle  precipitation  data.  The  boundary  migrates 
equatorward  with  increasing  magnetic  activity  [Dandekar,  1979],  Moreover, 
the  ionization  associated  with  the  diffuse  aurora  is  a  prominent  feature 
in  the  Chatanika  radar  data  [Vondrak  and  Baron,  1976], 

We,  therefore,  expect  the  diffuse  auroral  boundary  to  be  a  prominent 
feature  in  TEC  data  measured  by  using  polar  orbiting  satellites.  We  have 
found,  however,  that  the  F-region  ionization  associated  with  the  afore¬ 
mentioned  instability  [Vickrey,  1981]  dominates  the  TEC  boundary  structure. 
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B.  The  Occurrence  and  Location  of  TEC  Boundaries 

In  Figure  9,  we  show  two  nighttime  Wideband  satellite  TEC  records 
that  have  well-defined  TEC  boundaries.  The  data  are  presented  here  in 
units  of  relative  phase  change.  To  characterize  the  boundary,  the  mini¬ 
mum  and  maximum  locations  have  been  identified  (see  Figure  9)  together 
with  the  corresponding  change  in  TEC  units. 

As  discussed  in  Section  III-A,  the  presence  of  an  enhanced  F-layer 
complicates  the  interpretation  of  the  boundary.  As  a  guideline,  there¬ 
fore,  we  have  performed  model  computations  similar  to  those  used  in  Rino 
and  Owen  [1980b].  In  Figures  10,  representative  trough  and  diffuse  auru- 
ral  electron  density  profiles  are  shown  together  with  a  meridional  profile 
of  the  corresponding  density  contours  are  shown.  The  resulting  TEC  vari¬ 
ation  is  shown  in  Figure  11. 

The  results  of  a  similar  set  of  computations  for  an  enhanced  F  layer 
is  shown  in  Figures  12  and  13.  The  main  difference  is  that  the  F  region 
produces  a  much  larger  gradient.  The  location  of  the  boundary  relative 
to  the  station  is  also  important  as  can  be  seen  from  Figures  14  and  15. 

The  F-region  boundary  is  visible  from  approximately  70°  magnetic  latitude, 
whereas  the  E-region  boundary  is  only  visible  from  a  few  degrees  poleward 
of  the  station.  We  conclude  that  the  occurrence  of  a  well-defined  TEC 
boundary  is  heavily  weighted  toward  F-region  contributions,  particularly 
when  the  gradient  is  large  and  observed  poleward  of  the  station. 

With  this  in  mind,  we  consider  the  occurrence  of  TEC  boundaries  in 
the  nighttime  Wideband  passes  recorded  at  Poker  Flat  from  May  1976  through 
May  1979.  The  summary  data  are  shown  in  Figure  16.  A  TEC  boundary  is 
present  in  approximately  50  percent  of  the  passes,  irrespective  of  season 
or  solar-cycle  epoch. 

The  location  of  the  boundary  is  shown  in  Figure  17.  The  TEC  boundary 
moves  equatorward  with  increasing  activity,  but  generally  maintains  a 
constant  latitudinal  extent.  The  more  disturbed  passes  show  an  additional 
small  equatorward  migration  with  increasing  solar  activity.  The  three- 
year  period  from  1976  to  1979  is  generally  one  of  increasing  solar 
activity  as  the  1980  maximum  of  solar  cycle  21  was  approached. 
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MAGNETIC  LATITUDE  —  degrees 


FIGURE  9 


TYPICAL  NEAR-MERIDIONAL,  HIGH-LATITUDE 
TEC  VARIATIONS  SHOWING  BOUNDARY  FEATURE 
(ARROWS  INDICATE  MINIMUM  AND  MAXIMUM 
BOUNDARY  LOCATIONS) 


AURORAL  e-region 


TOTAL  ELECTRON  CONTENT  x  10  —  el/m' 


FIGURE  11  CALCULATED  MERIDIONAL  TEC  PROFILE  FOR  E-REGION  lONIZATlOf 
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FIGURE  12  UNIFORM  AURORAL  ENHANCED  F-REGION  IONIZATION  PROFILE 
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IGURE  14 


F  REGION 


effect  of  varying  enhanced  f-region 


RE  16  ANNUAL  AND  SEASONAL  OCCURRENCE  OF  TEC  BOUNDARY  FEATURES 


MAGNETIC  LATITUDE  —  degrees 


FIGURE  17  TEC  BOUNDARY  LOCATION  FOR  PERIODS  OF  LOW  AND  MODERATE 
LOCAL  MAGNETIC  ACTIVITY 


ATEC 


Structure  of  TEC  Boundary 


To  characterize  the  structure  of  the  TEC  boundary,  we  first  measured 
the  TEC  gradient.  In  Figure  18,  the  results  for  weakly  and  moderately 
disturbed  passes  are  shown.  The  gradients  evidently  depend  on  the  season, 
in  that  they  are  generally  larger  in  the  summer  months  than  in  the  winter 
months;  moreover,  the  TEC  gradient  generally  increases  with  solar  activity, 
as  does  the  occurrence  of  detectable  TEC  gradients. 


•4 -  1976  - ►  << -  1977  - ►  -  1978  - ^1979 


FIGURE  18  TEMPORAL  VARIATION  OF  TEC  GRADIENT  SHOWING  AN  INCREASE  WITH  INCREASING 
SOLAR-CYCLE  ACTIVITY 


41 


To  explore  the  seasonal  dependence  in  more  detail,  a  linear  least- 
squares  fit  to  the  TEC  gradient  as  a  function  of  the  corresponding  three- 
hour  college  Kindex  has  been  plotted  for  each  six-month  period.  The  re¬ 
sults  are  summarized  in  Figure  19.  The  gradient  decreases  rapidly  with 
increasing  K  during  the  fall  and  winter  months,  whereas  during  the  summer 
months,  it  increases. 

We  have  associated  the  steep  TEC  gradients  with  an  enhanced,  unstable 
F  layer.  To  show  this  directly,  we  have  scatter  plots  of  the  average  rms 
phase  measured  between  the  TEC  gradient  boundaries  against  the  measured 
gradient.  In  all  cases,  the  rms  phase  shows  a  positive  correlation. 

The  magnitude  of  the  slope  of  the  linear  least-squares  fit  for  each 
period  is  summarized  in  Table  1. 


1976  1977  1978  1979 


FIGURE  19  SLOPE  OF  LINEAR  LEAST-SQUARES  FIT  TO  TEC 
GRADIENT  vs  COLLEGE  K  INDEX 
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Table  1 


SLOPE  OF  LINEAR  LEAST  SQUARES  FIT  TO  rms  PHASE 
VERSUS  TEC  GRADIENT  BY  6 -MONTH  CALENDAR  PERIOD 


Year 

Period 

Slope 

1976 

1 

0.142 

2 

0.108 

1977 

1 

0.111 

2 

0.216 

1978 

1 

0.185 

0.154 

D.  Conclusions 

In  this  section  we  have  reviewed  the  occurrence  and  structure  of 
nighttime  TEC  boundaries  as  observed  by  the  Wideband  satellite  from  the 
Poker  Flat  Station  between  May  1976  and  June  1979.  Well-defined  TEC 
boundaries  occur  in  about  50  percent  of  the  passes  irrespective  of  the 
season  or  solar-cycle  epoch.  The  steepest  gradients  have  been  associated 
with  an  enhanced,  but  unstable  F  region.  The  occurrence  statistics  pri¬ 
marily  reflect  the  incidence  of  enhanced  F-region  "blobs"  migrating 
equatorward  of  70°  magnetic  latitude. 

The  TEC  boundary  structure  strongly  depends  on  magnetic  activity 
and  more  weakly  on  the  solar-cycle  epoch.  The  TEC  gradients  tend  to  be 
positively  correlated  with  magnetic  activity  during  the  summer  months, 
but  negatively  correlated  during  the  winter  months.  In  all  cases  the 
phase-scintillation  level  measured  within  the  TEC  boundary  region  is  posi¬ 
tively  correlated  with  the  TEC  gradient. 
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V  EQUATORIAL  TEC 


A.  Introduction 

The  Wideband  coherent  beacon  observations  provide  not  only  a  charac¬ 
terization  of  ionospheric  phase  scintillation,  but  also  a  measure  of  the 
large-scale  (>  100  km)  variations  that  are  directly  related  to  ionospheric 
TEC.  In  this  section  we  consider  the  TEC  measurements  made  at  Kwajalein 
Atoll  and  Ancon,  Peru,  both  equatorial  stations,  between  May  1976  and 
August  1979.  We  first  consider  the  nature  of  the  TEC  measurement,  and 
the  difficulties  encountered  in  applying  the  most  obvious  methods  for 
interpreting  the  data.  A  brief  review  of  equatorial  electrodynamics  and 
studies  of  three  aspects  of  our  TEC  observations  follow:  (1)  the  dramatic 
sunspot  number  dependence  of  the  midnight  equatorial  TEC,  (2)  the  seasonal 
behavior  of  the  latitudinal  distribution  of  TEC,  and  (3)  the  association 
between  TEC  and  the  occurrence  of  scintillation-producing  irregularities. 


Analysis  Method 


1.  Measurements 


Each  of  the  dispersive  phase  measurements  available  at  nine 
frequencies  from  the  Wideband  system  [Fremouw  et  al.,  1978]  is  directly 
related  to  the  propagation-path-integrated  electron  density 
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where  b  is  a  constant,  f  is  the  measurement  frequency,  and  s  indicates 
the  propagation  path  to  the  satellite.  In  practice,  the  dispersive  phase 
is  continuously  reconstructed  from  high  sample  rate  quadrature  components, 
and  thus  its  resolution  is  limited  only  by  the  receiver,  i.e.,  changes  in 
the  order  of  degrees  can  be  observed.  However,  because  of  2ir  ambiguities, 
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the  reconstruction  has  an  unknown  starting  value,  which  for  absolute 
measure  of  TEC  must  be  established  by  some  other  means. 

In  theory,  the  UHF  comb  of  seven  frequencies  on  Wideband  can 
also  provide  the  necessary  absolute,  but  lower  resolution,  measure  of 
TEC.  The  method  uses  the  second-difference  of  phases  for  a  number  of 
closely  spaced  spectral  lines  [Burns  and  Fremouw,  1970] 


f  2 

A,cp  b  — T 
£  f 


s 


Nds 


(12) 


where  f  is  the  frequency  separation  in  the  spectral  comb  (i.e.,  11.4729 
MHz  for  Wideband). 

During  the  data  collection  portion  of  the  experiment  much  effort 
was  expended  in  the  application  of  Eq.  (12).  This  was  necessary  because 
as  an  absolute  dispersion  measure,  is  sensitive  not  only  to  ionospheric 
dispersion,  but  any  dispersion  in  the  spacecraft  transmission  or  ground 
receiver  systems.  The  methods  of  prelaunch  calibration  and  post-launch 
comparisons  with  other  measures  of  TEC  in  the  auroral  zone  are  detailed 
in  Fremouw  et  al.  [1977],  Similar  calibrations  were  made  at  Kwajalein 
using  ALTAIR  radar  Lincoln  Lab  Sphere  tracks,  and  at  all  three  ground 
stations  using  standardized  antennas  [Cousins,  1979],  The  resulting  sets 
of  system  dispersion  constants,  estimated  for  each  station,  were  applied 
during  the  processing  of  each  satellite  pass  to  supply  a  starting  value 
to  the  continuous  reconstruction  of  high  resolution  dispersive  phase. 

During  the  subsequent  compilation  of  the  TEC  data,  on  average, 
the  TEC  values  seem  correct;  however,  pass-to-pass  variations  are  often 
much  too  large  to  be  attributed  to  ionospheric  changes.  Following  recent 
detailed  reanalysis  of  both  the  method  and  the  data,  we  have  concluded 
that  such  variations  are  unavoidable,  and  that  routine  calibrations  using 
A 2<P  are  not  practical. 

The  spacecraft  antenna  is  highly  frequency-dispersive,  because 
of  the  limited  size  of  the  ground  plane  and  multipath  effects.  The  pat¬ 
terns,  taken  before  launch  for  later  calibration  use,  are  highly  structured. 
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and  generally  subject  to  question,  since  other  aspects  of  the  patterns 
were  observed  to  change  after  launch.  The  UHF  antenna  was  least  dispersive 
at  low  angles,  and,  accordingly,  the  application  of  the  measurement 
was  made  early  in  each  pass.  A  second  conflicting  problem  is  introduced, 
in  that,  at  these  low  angles,  the  phase  scintillation  tends  to  maximize. 

The  combined  effect  of  these  problems,  i.e.,  noisy  data  and  dispersive 
structure  in  the  spacecraft  antenna,  contributes  errors  on  the  order  of 
tens  of  degrees  in  measured  ^  cp.  This  error  is  comparable  to  the  iono¬ 
spheric  contribution  at  the  most  closely  spaced  lines.  Only  in  a  few 
cases  has  it  been  possible  to  eliminate  obviously  noisy  data  and  improve 
the  TEC  estimate. 

Thus,  an  absolute  calibration  of  the  bulk  of  the  Wideband  TEC 
data  is  not  feasible.  Fortunately,  this  is  not  a  major  problem  for  a 
morphological  study  such  as  this,  because  the  shape  of  the  latitudinal 
TEC  distribution  reflects  the  solar  cycle  and  seasonal  dependences.  When 
the  term  TEC  is  used  to  describe  our  data  in  the  following  sections,  we 
are  actually  referring  to  changes  in  TEC,  as  derived  from  the  dispersive 
phase  alone. 

2.  Geometrical  Corrections 

The  high  inclination  of  the  Wideband  satellite  orbit  provides 
a  "snapshot"  of  the  ionosphere  over  about  20°  of  latitude  in  about  15  min. 
The  propagation  geometry,  of  course,  changes  significantly  from  overhead 
to  the  horizons,  and  interpretation  of  the  TEC  data  in  terms  of  the  lati¬ 
tudinal  distribution  of  ionization  requires  removal  of  that  propagation 
geometry.  The  simplest  and  most  widely  used  approach  is  to  assume  a  slab 
of  constant  ionizat  on  at  a  particular  height.  In  that  case,  the  TEC 

variation  can  be  described  simply  as  N_  =  N  sec0,  where  9  is  the  zenith 

1  o 

angle  of  the  propagation  path  at  the  chosen  layer  altitude,  and  Nq  is  the 
overhead  TEC  value. 

We  have  found  that  such  a  constant  altitude  slab  model  is  not 
consistent  with  the  Wideband  equatorial  data,  and  can  be  demonstrated 
using  some  of  the  early  Wideband  data  collected  in  1976  at  Ancon.  During 
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this  low  sunspot  phase  of  the  solar  cycle  (SSN  ~  15),  the  evening  F-region 
vertical  drifts  are  small,  and,  therefore,  the  equatorial  anomaly  in  TEC 
is  absent  at  the  time  of  Wideband  observation,  which  is  around  local  mid¬ 
night  (see  Section  C,  below).  During  this  period,  the  shape  of  the  un¬ 
corrected  TEC  records  are  nearly  constant  (for  the  same  propagation  geom¬ 
etry)  for  weeks  at  a  time.  Therefore,  application  of  an  appropriate 
geometry  correction  scheme  should  produce  similar  latitudinal  variations 
in  TEC  for  both  high-  and  low-elevation  passes. 

The  application  of  the  constant  height  slab  for  these  data  is 
shown  in  Figure  20  for  two  heights  and  several  passes  with  different 
propagation  geometries.  These  have  been  grouped  into  two  sets:  those 
with  maximum  elevation  angles  of  30°  to  60°  for  which  the  sec  9  factor 
changes  from  —4  to  ~  2  from  the  horizons  to  closest  approach,  and  those 
with  maximum  elevation  angle  of  60°  to  90°  for  which  sec  0  changes  from 
~  four  to  1.  If  the  constant  height  slab  model  is  appropriate,  the  shapes 
of  these  geometry  corrected  distributions  should  be  the  same  (although 
not  necessarily  flat).  This  is  clearly  not  the  case  for  a  layer  between 
350  km  and  700  km  in  Figure  20(a)  and  20(b);  the  slab  model  produces  a 
rapid  change  in  TEC  around  closest  approach,  which  differs  for  the  two 
pass  geometries. 

A  model  that  produces  a  more  consistent  geometrical  correc¬ 
tion,  is  one  for  which  the  F-layer  height  follows  a  particular  magnetic 
field  line.  Tsunoda  [1980]  has  presented  incoherent-scatter  data  from 
Kwajalein  that  clearly  shows  this  behavior  for  the  nighttime  F  layer.  In 
Figure  20(c),  we  show  the  same  low-high-elevation  pass  sets  from  Ancon, 
now  corrected  for  a  field-aligned  layer.  The  height  of  the  layer  ranges 
from  about  700  km  over  the  dip  equator  to  about  300  km  at  10°  dip  latitude. 
The  shapes  of  the  low-  and  high-elevation  latitudinal  distributions  are 
much  more  consistent  than  for  the  constant  height  layer  in  Figure  20(a) 
or  20(b). 

Having  established  a  geometrically  consistent  model  with  which 
to  interpret  the  data,  we  now  consider  the  latitudinal  variation  in  ef¬ 
fective  layer  thickness  that  produces  the  curvature  in  the  Figure  20(c) 
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data.  From  that  figure,  clearly,  the  curvature  cannot  be  removed  by 
moving  the  layer  up  or  down  (i.e.,  by  the  sec  0  dependence  on  layer  alti¬ 
tude).  This  idea  is  further  supported  by  the  ionosonde  data  from  Kwaja- 
lein,  which,  despite  an  increase  in  fQF2  from  1977  to  1979,  showed  little 
long-term  change  in  the  height  of  the  F-region  peak. 

The  most  direct  evidence  that  the  curvature  in  Figure  20  is  due 
to  a  change  of  layer  thickness  with  latitude,  however,  is  that  both  the 
Ancon  and  Kwajalein  data  show  similar  patterns  during  this  period  of  low 
solar  activity.  Because  only  a  few  days  of  data  from  Kwajalein  are  avail¬ 
able  during  the  September-to-October  period,  we  illustrate  the  similarity 
between  stations  using  the  November-to-December  1976  data.  As  Figure  21 
shows,  the  ionization  distribution  has  become  slightly  asymmetrical  during 
this  season,  because  of  the  equatorial  anomaly  to  the  north  (see  Section  C, 
below).  The  shapes  of  the  geometrically  corrected  TEC  distributions  with 
latitude  are  very  similar  despite  the  large  longitudinal  separation  of 
the  two  stations.  We  feel  that  this  justifies  the  use  of  the  Ancon 
September-to-October  distribution  as  our  low  sunspot  number,  no-anomaly 
baseline  for  the  latitudinal  distribution  corrections  at  Ancon  and 
Kwajalein. 

Variation  in  F-layer  thickness  with  latitude  is  also  consistent 
with  the  single  data  set  presented  in  Tsunoda  [1980],  in  which  he  observed 
much  narrower  F-layer  thickness  with  roughly  the  same  peak  density  south 
of  Kwajalein.  (His  observations  to  the  north  of  Kwajalein  cannot  be  com¬ 
pared  to  our  1976  data  because  of  higher  sunspot-number  conditions.) 

For  study  of  the  long-term,  large-scale  nighttime  TEC  behavior, 
then,  we  postulate  the  following  low  sunspot  background  model:  an  F  layer 
whose  centroid  is  field  aligned,  ranging  in  altitude  from  700  km  over  the 
dip  equator  to  about  350  km  at  10°  dip  latitude,  and  an  effective  thick¬ 
ness  at  the  magnetic  equator  roughly  double  that  at  ±10°  dip  latitude. 

In  what  follows,  we  consider  the  long-term  deviations  from  this  baseline 
model  because  of  sunspot  number,  season,  and  other  causes.  Before  the 
data  are  presented,  however,  a  brief  description  of  the  general  electro¬ 
dynamic  behavior  of  the  equatorial  ionosphere  as  it  is  currently  under¬ 
stood  is  worthwhile. 
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FIGURE  21  GEOMETRICALLY  CORRECTED  TEC  OBSERVATIONS  AT 
(a)  ANCON,  AND  (b)  KWAJALEIN  FOR  NOVEMBER 
THROUGH  DECEMBER  1976  LOW  SUNSPOT  CONDITIONS 


C.  Dynamics  of  the  Equatorial  F-Region 

The  behavior  of  the  equatorial  F  region  is  dominated  by  electric 
fields.  Foremost  of  these  is  the  E-region  polarization  field  created  by 
tidally  controlled  winds  blowing  across  the  geomagnetic  field.  Fields 
generated  away  from  the  magnetic  equator  map  directly  into  the  F-region 
along  the  highly  conducting  flux  tubes.  This  E-region  dynamo  drives  the 
F-region  ionization  continuously  upward  by  E  X  B  forces  during  the  day. 
From  these  high  altitudes,  pressure  and  gravity  force  the  ionization  back 
down  flux  tubes  to  higher  and  lower  latitudes.  This  rise  and  subsequent 
downward  diffusion  produces  a  latitudinal  variation  in  density  that  peaks 
to  the  north  and  south  relative  to  a  trough  over  the  equator:  the  "Apple' 
ton  anomaly." 

Many  years  of  ionosonde  measurements  have  defined  the  overall  depend, 
ence  of  the  anomaly  on  sunspot  number,  season,  and  time  of  day.  During 
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Che  day,  at  the  equinox  and  ac  sunspot  minimum,  the  crests  of  the  anomaly 
occur  at  ~17°  dip  latitude  with  a  critical  frequency  roughly  double  that 
at  the  equator  [Rastogi,  1959].  A  similar  configuration  has  been  reported 
by  Lyon  and  Thomas  [1963]  for  the  same  time  and  season,  but  sunspot  maxi¬ 
mum.  The  sunspot  number  lias  a  strong  effect  on  the  duration  of  the  ano¬ 
maly  into  the  evening  hours,  which  are  of  interest  here.  Rastogi  _1959_ 
shows  that  for  sunspot  minimum  that  the  anomaly  disappears  by  2100  LT. 

On  the  other  hand,  Lyon  and  Thomas  [1963_  show  that  the  anomaly  becomes 
more  distinct  in  the  premidnight  hours  during  sunspot  maximum.  At  the 
solstices,  north-south  asymmetries  in  f^l’2  typically  occur  ]Lyon  and 
Thomas,  1903]  with  a  reversaL  of  that  asymmetry  from  day  to  night. 

Modeling  of  these  seasonal  and  diurnal  variation  of  the  anomaly  as 
driven  by  the  E-region  dynamo  have  been  quite  successful.  !:or  example, 
Anderson  ]  1 9 7 3 _  solves  the  time-varying  plasma  continuity  equation  in¬ 
cluding  the  effects  of  magnetic  field,  neutral  wind,  and  E-region  dynamo 
fields.  Using  that  model,  he  duplicates  the  experimentally  observed  sua- 
tial  and  diurnal  patterns  of  the  anomaly. 

However,  certain  key  features  of  the  <.  zoning  i-'- layer  behavior  cannot 
be  explained  by  the  E-region  dynamo  alcu.e.  The  rapid  rise  of  the  F-region 
at  sunset,  indicating  a  strongly  enhanced  eastward  electric  field,  is  a 
case  in  point.  Rishbeth  [1971]  has  explained  this  behavior,  which  is  re¬ 
ferred  to  as  the  prereversal  enhancement,  using  an  F-region  dynamo  driven 
by  thermospheric  winds.  During  the  day,  these  fields  are  shorted  by  a 
highly  conducting  E-region,  but  near  sunset,  while  high  residual  F-region 
ionization  remains  and  strong  neutral  winds  are  present,  F-region  genera¬ 
ted  fields  become  significant.  Heelis  et  a 1 ,  [l974j  have  modeled  this 
process  with  good  agreement  with  observation  at  Jicamarca.  Fejer  et  al. 
[1979]  have  compiled  a  large  volume  of  vertical-drift  data  from  Jicamarca 
taken  over  several  years.  Their  evening  and  nighttime  data,  which  are 
of  interest  to  this  study,  are  dominated  by  the  enhancement  in  vertical 
drift,  which  depends  highly  on  both  sunspot  number  and  season.  Further¬ 
more,  the  time  of  reversal  from  upward  to  downward  driftf  in  the  evening 
shows  similar  dependences. 
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In  addition  to  the  electric  fields  produced  by  F-region  zonal  winds, 
the  meridional  component  of  the  wind  has  a  direct  effect  on  the  latitudi¬ 
nal  distribution  of  ionization  [Rishbeth,  1977].  The  winds  flow  from 
]  high  to  low  pressure  regions  which  have  been  observed  to  be  highly  asymme¬ 

trical  across  the  equator  [Hedin  et  al..  1974],  In  general  the  flows  are 
from  the  summer  to  the  winter  hemisphere,  although  the  diurnal  and  short¬ 
term  (few  days)  perturbations  to  this  pattern  can  be  significant. 

The  effect  of  the  transequatorial  wind  transfers  plasma  from  the  up¬ 
wind  side  of  the  equator  to  the  downwind  side,  therefore  producing  sig¬ 
nificant  asymmetries  in  the  equatorial  anomaly.  Anderson  and  Matsushita 
[1974]  have  estimated  the  order  of  magnitude  of  this  variation  with  season. 
At  times,  diurnal  variations  in  neutral  pressure  tend  to  produce  diverging 
(day)  or  converging  (night)  rather  than  transequatorial  winds.  The  ano¬ 
maly  under  such  conditions  has  been  modeled  by  Rush  [1972], 

D.  Wideband  A  TEC  Results 

1.  Overall  Behavior  and  Sunspot  Dependence 

In  Figure  22,  three  years  of  TEC  data  from  Ancon  are  shown,  the 
1976  data  at  the  top  and  1979  data  at  the  bottom.  Each  frame  corresponds 
to  two  months  of  observations,  from  which  the  propagation  geometry  and 
the  low-sunspot  latitudinal  distribution  of  Figure  20c  have  been  removed. 
Furthermore,  because  we  have  no  absolute  density  calibration  for  each 
pass,  all  have  been  assigned  an  appropriate  constant  value  of  TEC  at  the 
magnetic  equator,  where  over  the  long  term,  the  nighttime  overhead  TEC 
should  least  depend  on  electrodynamic  processes. 

Our  measurements  extend  to  about  ±11°  dip  latitude,  which  is 
well  into  the  region  in  which  the  equatorial  anomaly  crests  of  fQF2  nor¬ 
mally  occur.  Because  TEC  is  an  integrated  measure  of  the  F-layer,  it  is, 
in  some  ways,  a  more  complete  indicator  of  ionization  deposition  in  the 

anomaly  crests  than  f  F2  measurements.  However,  if  the  thickness  of  the 

o 

F-layer  under  anomaly  conditions  changes  rapidly  with  latitude  [as  seen 
by  Vila,  1966,  for  daylight  hours],  the  TEC  distribution  with  latitude 
may  differ  from  that  implied  from  fQF2  studies. 
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As  Figure  22  shows,  the  midnight  ionosphere  at  Ancon  changed 
very  little  during  the  low  sunspot  conditions  of  July  through  October 
1976.  Corresponding  to  the  upswing  of  the  solar  cycle  from  late  1976 
through  the  end  of  our  observations,  we  see  a  dramatic  increase  in  TEC 
magnitudes  in  the  anomaly  crest  region. 

The  general  shapes  of  the  ionization  distributions  in  Figure  22 
show  that  during  the  spring  equinox,  the  Ancon  observations  extend  through 
the  southern  crest  of  the  anomaly.  This  is  most  conspicuous  in  the  1977 
and  1978  data,  and,  although  a  similar  crest  is  not  clear  in  the  1979 
data,  we  can  make  a  conservative  estimate  of  the  increase  of  TEC  in  the 
anomaly  crests  with  increasing  sunspot  number.  The  distributions  imply 
that  under  equinox  conditions  the  anomaly  TEC  from  1977  and  1978  changes 
~60%  and  at  least  another  60%  from  1978  to  1979,  i.e.,  more  than  a  doub¬ 
ling  in  TEC  over  the  first  half  of  the  solar  cycle  upswing. 

As  the  sunspot  number  increases,  so  does  the  magnitude  of  the 
TEC  in  the  anomaly  crest  region,  and  the  overall  TEC  variability.  Much 
of  the  variability  is  caused  by  systematic  solar,  seasonal,  and  diurnal 
changes  in  bulk  ionization  transport.  During  the  scintillation  season 
our  data  also  includes  the  signatures  of  specific  ionospheric  features, 
i.e.,  Rayleigh-Taylor  bubbles,  the  occurrence  of  which  also  increases 
with  solar  activity.  The  high  inclination  of  the  Wideband  orbit  brings 
the  satellite  nearly  along  the  magnetic  meridian  and  if  a  bubble  aligned 
with  a  particular  flux  tube  is  intercepted,  it  may  remain  in  the  propaga¬ 
tion  path  for  several  minutes  of  the  entire  ~  15-min  north-south  scan. 

The  distinct  signatures  of  such  bubbles  can  be  seen  in  some  of  the  early 
1978  and  1979  data  in  Figure  22. 

The  Ancon  data  can  be  directly  compared  to  the  vertical  drift 
data  collected  at  the  Jicamarca  radar  [Fejer  et  al.t  1979].  By  comparing 
Figure  22  with  those  data,  we  can  see  that  the  dramatic  change  in  anomaly 
TEC  with  sunspot  number  is  closely  associated  with  the  F-region  dynamo, 
or  more  specifically  with  the  enhanced  evening  rise  of  the  ionosphere. 
Fejer  et  al.  _1979j  have  shown  that  the  magnitude  of  that  rise,  which  may 
exceed  30  m/s,  depends  highly  on  solar  cycle.  At  sunspot  minimum,  the 
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prereversal  enhancement  is  observed  only  at  the  equinoxes,  and  even  then 
is  weak.  This  is  consistent  with  our  1976  data,  although  we  see  no  evi¬ 
dence  of  any  anomaly  even  in  the  Fall  1976  equinox,  which  probably  im¬ 
plies  that  excess  ionization  transported  by  the  weak  prereversal  enhance¬ 
ment  had  decayed  by  local  midnight. 

Fe/jer  et  al.  [1979]  also  find  that  during  higher  solar  activity 
conditions  the  prereversal  enhancement  is  strongest  in  March  and  September; 
in  1977  we  see  an  anomaly  peak  for  the  same  months.  The  enhancement  is 
weakest  from  May  to  August,  but  does  have  a  longer  duration,  and  continues 
until  ~  2100  LT,  some  two  hours  later  than  in  the  winter.  We  observe  much 
less  anomaly  from  May  to  August  1977  than  in  the  equinoxes,  which  implies 
that  primarily,  the  magnitude  of  the  enhancement,  and  not  its  duration, 
is  important  to  the  formation  of  the  nighttime  anomaly.  This  idea  is  not 
too  surprising  because  the  lower  ionosphere  densities,  from  which  the 
anomaly  is  formed,  are  still  high  at  E-region  sunset,  then  rapidly  decrease. 

Turning  to  the  Kwajalein  data  in  Figure  23,  we  see  a  somewhat 
different  latitudinal  distribution  of  TEC.  Because  the  observing  station 
at  Kwajalein  is  at  ~  5°  dip  latitude,  we  expected  to  see  more  of  the 
northern  portion  of  the  anomaly.  Figure  23  shows  that  this  is,  indeed, 
the  case  and  in  fact,  the  full  crest  of  the  TEC  anomaly  is  seen  consist¬ 
ently  at  12°  north  dip  latitude. 

A  direct  comparison  can  be  made  between  the  Kwajalein  and  Ancon 
data  for  the  periods  of  simultaneous  observations.  For  the  low  sunspot 
number  conditions  of  November  to  December  1976,  the  TEC  distributions  at 
the  two  stations  are  very  similar,  reassuring  us  that  the  model  layer  re¬ 
moval  discussed  in  Section  V-B  is  appropriate.  The  January-to-February 
1977  data  are  more  variable  at  Ancon  than  at  Kwajalein,  but  the  median 
TEC  distributions  for  this  and  the  remaining  five  periods  of  simultaneous 
observations  match  closely.  During  this  period  of  moderate  solar  activity 
then,  the  large-scale  electrodynamics  processes  producing  the  anomaly 
must  be  the  same  in  the  American  and  Pacific  sectors,  despite  their  large 
longitudinal  separation. 
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FIGURE  23  LATITUDINA 


Unlike  the  Ancon  data,  which  spans  September  through  April  and 
is  concentrated  near  the  equinoxes,  the  Kwajalein  observations  were  made 
primarily  during  May  through  August.  These  are  the  months  during  which 
the  prereversal  enhancement  in  vertical  drift  is  minimum.  As  we  did  at 
Ancon,  we  see  direct  evidence  of  that  seasonal  change  in  the  Kwajalein 
data.  Between  May  and  August  of  both  1977  and  1978,  we  observe  a  constant 
TEC  value  at  the  summertime  northern  anomaly  crest  (and  actually  a  de¬ 
crease  in  1978),  when  compared  to  the  equinox  value.  This  seasonal  varia¬ 
tion  is  further  verification  of  F-region  dynamo  control.  The  relatively 
constant  magnitude  of  TEC  in  the  northern  anomaly  crest  between  the  summers 
of  1978  and  1979,  despite  a  doubling  of  the  sunspot  number,  is  less  . asy 
to  explain.  This  compares  with  the  equinox  TEC  increase  of  >  60  percent 
between  the  same  years  as  observed  at  Ancon.  Unfortunately,  we  do  not 
have  enough  equinox  data  from  Kwajalein  to  guarantee  that  this  is  not  a 
longitudinal  effect,  although  no  evidence  exists  in  the  1977  and  1978 
eqj inox  data  that  shows  such  a  difference.  More  likely,  the  strengthen¬ 
ing  of  th.  F-region  dynamo  with  solar  activity  saturates  at  moderate  sun¬ 
spot  numbers  between  May  and  August.  No  such  saturation  can  be  implied 
from  the  Ancon  data  for  the  November-through- February  solstice  period. 

We  next  consider  the  seasonal  changes  of  the  TEC  distributions 
in  more  detail. 

2.  Seasonal  Dependence 

In  Figures  22  and  23,  the  envelopes  of  the  TEC  distributions 
indicate  a  large  degree  of  variability  within  each  two-month  period. 

Some  of  the  variability  is  caused  by  the  intersection  with  discrete, 
medium-scale  structures  (i.e.,  bubbles)  as  noted  above.  Much  of  it,  how¬ 
ever,  is  caused  by  systematic  trends  with  periods  ranging  from  a  few  days 
to  a  few  weeks.  In  this  section,  we  consider  the  data  in  terms  of  those 
trends,  and  speculate  on  their  cause. 

Starting  again  with  the  Ancon  data,  Figure  24  shows  contours  of 
constant  TEC  versus  latitude  and  day  number  for  the  first  half  of  1977, 
1978,  and  1979.  Each  contour  level  has  been  assigned  a  TEC  value  based 


FIGURE  24  TEC  VARIATION  WITH  LATITUDE  AND  SEASON  AT  ANCON,  JANUARY  THROUGH 
JUNE,  1977  THROUGH  1979 
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on  a  typical  overhead  value  of  2.2  X  10  el/m  (150  rad  dispersive  phase 
at  1239  MHz).  Some  weeks  include  more  passes  than  others,  which,  of 
course,  contributes  to  the  structure  in  the  contours.  We  are  concerned 
here  with  the  overall  behavior;  in  the  next  section,  when  we  compare  the 
TEC  and  scintillation  occurrence,  more  detail  will  be  considered. 

In  Figure  24,  the  increase  in  TEC  with  increasing  sunspot  ac¬ 
tivity  is  clear,  as  in  Figure  22.  For  the  1977  data,  the  weakening  of 
the  anomaly  in  May  and  June  is  also  conspicuous  despite  the  continuous 
increase  in  sunspot  number  during  this  period.  However,  Figure  24  addi¬ 
tionally  shows  a  distinct  and  systematic  northward  shift  in  the  anomaly 
between  mid-March  and  mid-April  for  all  three  years.  In  1977,  the  anomaly 
was  to  the  north  and  was  strongest  near  Day  40;  in  1978,  the  occurrence 
peak  is  more  diffuse,  showing  subpeaks  at  about  Days  50  and  70;  in  1979, 
the  northern  anomaly  peaks  at  about  Day  50.  For  all  three  years,  the 
distribution  pattern  shifts,  and  the  northern  anomaly  disappears  and  is 
replaced  by  the  southern  peak  on  about  Day  94. 

In  Figure  25,  the  second  half  of  1976,  1977,  and  1978  are 
shown.  For  the  latter  two  years,  the  indication  is  that  the  northward 
shift  at  the  spring  equinox  is  reversed  in  the  fall.  Unlike  the  spring 
data,  however,  the  time  of  the  fall  southward  shift  appears  to  vary.  In 
1976,  no  real  anomaly  appears  to  the  south,  but  appears  weakly  to  the 
north  centered  on  Day  330.  In  1977,  the  anomaly,  which  appears  to  the 
south  on  about  Day  280,  shifts  to  the  north  by  Day  310.  In  1979,  a  strong 
southern  anomaly,  also  centered  on  Day  280,  is  replaced  by  a  weaker  ano¬ 
maly  to  the  north  some  20  to  30  days  later. 

Unfortunately,  except  for  the  low  sunspot  conditions  of  1976 
and  1977,  little  can  be  learned  from  this  data  set  about  the  spatial/ 
seasonal  behavior  between  late  December  and  February.  All  that  can  be 
implied  is  that  the  anomaly  weakens,  as  the  Jicamarca  vertical  drift  data 
predict . 

To  compare  the  seasonal  behavior  of  the  Ancon  and  Kwajalein 
data,  the  simultaneous  observation  times  during  1977  are  plotted  in 
Figures  26  and  27,  for  January  through  June  and  July  through  December, 
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FIGURE  25  TEC  VARIATION  WITH  LATITUDE  AND  SEASON  AT  ANCON,  JULY  THROUGH 
DECEMBER,  1976  THROUGH  1978 
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FIGURE  26  COMPARISON  OF  TEC  VARIATION  WITH  LATITUDE  AND  SEASON  AT  ANCON  AND  KWAJALEIN,  JANUARY  THROUGH 
JUNE  1977 


respectively.  As  can  be  seen,  the  anomaly  5'iifts  at  the  equinoxes  are 

not  as  clear  in  the  Kwajalein  data,  although  the  trends  are  the  same  as 
at  Ancon;  in  the  spring,  the  TEC  distribution  generally  shifts  northward 
by  about  2°,  very  much  like  that  shown  by  the  distribution  at  Ancon.  The 
southern  anomaly  appears  at  both  stations  on  or  about  Day  100.  In  the 
fall,  the  southward  shift  is  subtle,  but  as  might  be  expected,  the  north¬ 
ern  anomaly  appears  sooner  at  Kwajalein  than  at  Ancon. 

Similar  behavior  of  the  equinox  anomaly  can  be  seen  in  the 
spring  1978  data,  which  are  compared  in  Figure  28.  At  Ancon,  the  south¬ 
ern  anomaly  shifts  southward  after  a  peak  near  Day  96;  a  similar  shift 
southward  can  be  seen  in  the  Kwajalein  data.  During  the  shift,  the 
northern  crest,  of  which  there  is  some  evidence  at  Ancon,  also  progresses 
southward . 

The  biggest  difference  between  the  simultaneous  Kwajalein  and 
Ancon  data  sets  is  the  behavior  during  the  llay-through-August  period.  In 
Figures  22  and  23  the  TEC  distributions  generally  agreed  for  this  period. 
From  Figures  26  and  27,  we  see  now  that  a  match  occurs  only  because  the 
Ancon  measurements  are  concentrated  before  Day  150  and  after  Day  220, 
periods  when  the  contours  roughly  match.  Although  the  total  disappearance 
of  the  anomaly  between  Days  200  and  210  at  Ancon  cannot  be  confirmed, 
the  evidence  is  sufficient,  in  particular  the  anomaly  enhancements  at 
Kwajalein  near  Days  120,  145,  and  182,  that  the  anomaly  behavior  at  the 
two  stations  is  different.  This  suggests  that  a  longitudinal  difference 
exists  in  the  processes  creating  the  anomaly  under  summer  solstice  condi¬ 
tions  at  Kwajalein.  Unfortunately,  no  other  coincident  observations  be¬ 
tween  Kwajalein  and  Ancon  at  this  season  verify  or  refute  this  difference. 
Figure  26  shows  no  evidence  that  a  similar  difference  exists  during  the 
winter  solstice  at  Kwajalein  in  early  1977. 

From  the  accumulated  Kwajalein  data  the  northern  hemisphere 
summer  patterns  in  TEC  can  be  established  and  may  shed  some  light  on  the 
Kwajalein-Ancon  differences  noted  above.  The  contours  are  shown  in 
Figure  29  for  1977,  1978,  and  1979.  The  TEC  magnitudes  generally  in¬ 
creased  with  sunspot  numbers  from  1977  to  1978,  and  additionally,  but 


FIGURE  29  TEC  VARIATION  WITH  LATITUDE  AND  SEASON  AT  KWAJALEIN,  MAY  THROUGH  SEPTEMBER,  1977  THROUGH  1979 


with  a  lesser  degree,  increased  from  1978  to  1979,  consistent  with  Fig¬ 
ure  23.  No  patterns  are  obvious  in  the  anomaly  behavior  during  June  and 
July,  with  major  rearrangements  occurring  over  a  period  of  several  days. 
Correlation  of  particular  features  from  year  to  year,  such  as  that  at 
Day  180  are  most  likely  coincidental. 

The  year-to-year  variability  of  the  summer  data  and  the  differ¬ 
ences  between  Ancon  and  Kwajalein  during  this  season  both  point  to  local 
control  of  the  anomaly  pattern. 

E.  The  Association  Between  TEC  and  Scintillation 

Although  much  has  been  learned  about  the  equatorial  ionosphere  over 
the  past  few  years,  one  question  which  remains  is  the  seasonal/longitudinal 
control  of  scintillation.  Livingston  [1980]  has  shown  that  the  scintilla¬ 
tion  seasons  in  the  American  and  Pacific  sectors  are  roughly  eight  months 
long  and  six  months  out  of  phase.  Figure  30  illustrates  this  for  the 
1976-to-1977  intensity  scintillation  data.  The  season  centers  on  January 
in  the  American  (and  African)  sectors  and  July  in  the  Pacific.  Geographic 
control,  alone  can  be  discarded  as  a  cause  because  of  the  simultaneous 
scintillation  seasons  in  the  American  and  African  sectors.  A  combination 
of  geographic-geomagnetic  control  is  most  likely  the  controlling  factor, 
perhaps  through  the  E-region  dependent  process  modelled  by  Zalesak  et  al. 
[1981]. 

The  basic  tie  between  scintillation  and  TEC  is  that  the  scintillation 
depends  on  the  density  perturbation,  5n,  in  the  hundreds-of-meters  to  a 
few  kilometers  scale  regime;  the  higher  the  background  density,  the 
larger  6n  can  become.  With  this  in  mind,  we  have  compared  the  scintilla¬ 
tion  levels  with  the  TEC  observed  simultaneously  by  Wideband. 

The  severity  in  scintillation  generally  increases  at  both  stations 
between  1976  and  1977.  This  increase  can  be  attributed  directly  to  the 
solar-cycle-produced  increase  in  background  electron  density,  but  is  more 
likely  caused  by  a  combination  of  effects,  e.g.,  higher  density  plus  en¬ 
hanced  instability  growth  rate  from  F-layer  dynamic.  Any  direct  tie  in 
our  data  between  enhanced  TEC  and  scintillation  will  be  useful  in  clari¬ 
fying  the  specific  cause. 
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FIGURE  30  WIDEBAND  DATA  FROM  1977  SHOWING  SEASONAL  INTENSITY  SCINTILLATION 
OCCURRENCE  AT  la)  KWAJALEIN  AND  (b)  ANCON 


Detailed  TEC-scintillation  comparisons  have  been  made  for  campaign 
times  at  both  stations,  during  which  the  observations  were  more  or  less 
continuous.  To  characterize  the  scintillation,  we  have  chosen  the  stan¬ 
dard  deviation  of  the  phase  scintillation  at  413  MHz,  which,  in  this  case, 
defines  the  integrated  phase  spectral  energy  on  scales  shorter  than  ~  30  km. 
Like  the  TEC,  the  standard  deviation  of  phase  depends  on  the  effective 
thickness  of  the  ionospheric  layer  and,  thus,  must  be  corrected  for  geom¬ 
etry.  We  have  done  this  using  the  same  layer  height  as  used  in  the  TEC 
corrections. 

The  spring  equinox  Ancon  data  from  1977  are  shown  in  Figure  31a;  in 
Figure  31b,  the  corresponding  phase-scintillation  data  are  displayed, 
where  the  size  of  the  circle  indicates  the  magnitude  of  each  20-s  index. 
Values  less  than  0.2  rad  of  rms  phase  are  indicated  by  the  smallest  circles 
that  can  be  seen.  A  one-to-one  comparison  between  the  TEC  magnitudes  and 
the  scintillation  cannot  always  be  expected.  For  those  passes  that  inter¬ 
cept  a  depleted  flux  tube,  an  edge  of  which  may  be  structured,  the  oppo¬ 
site  can  be  true.  Thus,  we  should  look  for  trends  of  several  days  in 
the  occurrence  of  spatial  scintillation  to  match  the  same  trends  in  TEC. 

In  Figure  31,  there  is  clearly  no  evidence  of  a  consistent  corre¬ 
spondence  between  high  TEC  and  the  spatial/seasonal  scintillation  occur¬ 
rence.  The  anomaly  crests  early  to  the  north,  and  later  to  the  south  are 
not  associated  with  the  bulk  of  the  scintillation,  which  occurs  approxi¬ 
mately  between  Days  75  and  95.  The  scintillation  occurs  primarily  to 
the  south.  It  may  be  significant  that  this  three  weeks  of  peak  scintilla¬ 
tion  is  at  the  end  of  the  northward  shift  period  noted  in  Section  V-B. 

Interestingly,  the  same  pattern  occurs  in  the  1978  spring  data, 
shown  in  Figure  32.  Here  again,  association  between  TEC  and  scintilla¬ 
tion  is  not  clear.  In  fact,  during  the  span  from  Day  80  to  Day  100,  when 
the  anomaly  is  strongly  southward,  the  scintillation  is  conspicuously 
absent.  Instead,  as  in  Figure  31,  the  enhanced  scintillation  is  again 
primarily  southward,  and  occurs  during  the  northward  transition  of  the 
anomaly. 
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FIGURE  31  COMPARISON  OF  (a)  TEC  AND  (b)  UHF  PHASE  SCINTILLATION  WITH  LATITUDE  AND  TIME,  ANCON,  SPRING  1977 
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AND  TIME.  ANCON,  SPRING  1978 


We  have  also  considered  the  fall  period  in  1977,  when,  as  in  the 
spring,  the  Wideband  observations  were  concentrated.  These  data  are  shown 
in  Figure  33.  As  noted  in  Section  V-B,  the  fall  southward  shift  in  the 
anomaly  occurs  between  Days  280  and  320.  And,  as  Figure  33  shows,  the 
scintillation  occurrence  peaks  during  this  period,  rather  than  being 
associated  with  times  and  regions  of  enhanced  TEC.  These  consistent 
associations  between  scintillation  and  the  anomaly  shift  strongly  suggests 
the  influence  of  the  neutral  wind. 

The  Kwajalein  TEC-scintillation  data  for  1978  are  shown  in  Figure  34. 
First  note  the  general  seasonal  trend  in  average  scintillation.  The 
severity  of  the  peak  scintillation  in  April  is  comparable  to  that  in  July. 
The  average  scintillation  in  April  is  lower,  however,  because  of  the 
number  of  interspersed  quiet  passes.  The  overall  latitudinal  distribution 
of  scintillation  is  different  from  Ancon  in  that  the  disturbed  regions 
lie  to  the  north,  i.e.,  in  the  vicinity  of  the  anomaly. 

In  general,  the  l&rajalein  data  show  a  closer  association  between  en¬ 
hanced  TEC  and  scintillation  than  do  the  Ancon  data.  This  is  particularly 
true  during  the  period  after  Day  202  and  beyond.  The  exceptions  are 
notable,  however,  during  the  appearance  of  the  strongest  anomaly  between 
Day  100  and  Day  140,  and  the  moderate  anomaly  crests  between  Day  170  and 
Day  190. 

Figure  35  shows  a  more  disturbed,  but  shorter  span  for  the  1979 
Kwajalein  data.  The  patterns  for  the  1978  and  1979  scintillation  are 
not  consistent;  the  period  of  strongest  average  scintillation  is  earlier 
in  1979,  between  Day  175  and  Day  185.  The  correlation  between  scintilla¬ 
tion  and  TEC  in  1979  is  about  the  same  as  in  1978.  Both  Kwajalein  data 
sets  indicate  that  enhanced  TEC  regions  are  necessary  for,  but  do  not 
guarantee,  the  existence  of  scintillation  producing  irregularities. 

Significant  difference  exists,  then,  between  Ancon  and  Kwajalein  in 
terms  of  the  association  between  TEC  and  scintillation.  It  is  not  too 
surprising,  in  light  of  the  six-month  shift  in  scintillation  season  be¬ 
tween  the  American  and  Pacific  sectors,  that  different  mechanisms  may 
control  the  generation  of  scintillation  irregularities. 
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FIGURE  33  COMPARISON  OF  (a)  TEC  AND  (b)  UHF  PHASE  SCINTILLATION  WITH  LATITUDE  AND  TIME, 
ANCON,  AUTUMN  1977 


FIGURE  34  COMPARISON  OF  (a)  TEC  AND  (b)  UHF  PHASE  SCINTILLATION  WITH  LATITUDE  AND  TIME,  KWAJALEIN 
SUMMER  1978 


FIGURE  35  COMPARISON  OF  (a)  TEC  AND  <b>  UHF  PHASE  SCINTILLATION  WITH  LATITUDE 
AND  TIME,  KWAJALEIN,  SUMMER  1979 


F.  Conclusions 

Several  conclusions  can  be  drawn  from  our  equatorial  TEC  study. 

(1)  The  latitudinal  distribution  of  TEC  under  no-anomaly  con- 
disions  is  best  described  by  a  layer  whose  height  follows 
the  local  magnetic  field.  The  effective  vertical  TEC  of 
the  layer  at  ±10°  dip  latitude  is  roughly  half  that  at 
the  dip  equator,  implying  a  decreasing  layer  thickness  or 
peak  density,  or  both. 

(2)  The  occurrence  of  the  nighttime  equatorial  anomaly  is 
closely  tied  to  the  prereversal  enhancement  in  vertical 
drift  created  by  F-region  dynamo  fields.  The  anomaly  is 
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absent  during  low  solar  activity,  but  increases  rapidly 
with  increasing  sunspot  number,  as  does  the  prereversal 
enhancement.  Furthermore,  the  anomaly  is  strongest  at 
the  equinoxes,  and  weakest  during  the  northern  hemi¬ 
sphere,  again  in  accord  with  the  drift  enhancement.  This 
basic  electrodynamic  control  is  observed  in  both  the 
American  and  Pacific  sectors. 

(3)  The  northern  crest  of  the  anomaly,  as  observed  in  TEC, 
remains  within  a  few  degrees  of  12°  dip  latitude,  as  seen 
from  Kwajalein.  The  southern  crest  remains  generally  be¬ 
yond  our  observing  range  of  -11°  at  Ancon.  At  both  sta¬ 
tions,  a  distinct  northward  shift  in  the  anomaly  is  ob¬ 
served  following  the  spring  equinox.  A  similar,  but  less 
distinct,  southward  shift  accompanies  the  fall  equinox  at 
Ancon.  Such  shifts  in  the  anomaly  symmetry  are  a  result 
of  transequatorial  transport  of  ionization  by  F-region 
meridional  winds.  From  our  observations,  which  are 
limited  in  latitude,  it  cannot  be  established  whether 
this  shift  is  part  of  a  longer  term  pattern,  or  is  a 
distinct  equinoctal  event. 

(4)  While  the  seasonal  behavior  of  TEC  at  Ancon  and  Kwajalein 
is  similar  at  the  equinoxes,  it  differs  significantly 
during  the  northern  hemisphere  summer.  The  anomaly  disa¬ 
ppears  at  Ancon,  but  remains  distinct,  although  limited  in 
magnitude,  at  Kwajalein.  Also,  during  this  season,  the 
solar-activity  control  of  TEC,  observed  strongly  at  the 
equinoxes,  saturates  early  in  the  solar-cycle  upswing. 

(5)  The  association  between  the  occurrence  of  scintillation 
and  TEC  anomaly  is  not  conspicuous.  At  Ancon,  the  peak 
scintillation  is  apparently  associated  more  with  the 
equinox  anomaly  shift,  than  with  the  anomaly  magnitude. 

At  Kwajalein,  where  the  scintillation  season  is  shifted 
six  months  from  that  at  Ancon,  the  correlation  between 
TEC  magnitude  and  scintillation  is  better.  Enhanced  TEC 
is  generally  necessary  for,  but  does  not  guarantee,  en¬ 
hanced  scintillation. 

Overall,  then,  similarities  in  the  TEC  behavior  at  Ancon  and  Kwaja¬ 
lein  verify  the  basic  electric  field  dependence  of  bulk  motion  in  the 
equatorial  F  region.  On  the  other  hand,  differences,  such  as  the  sol¬ 
stice  behavior  of  the  anomaly,  and  the  scint  llation  behavior  at  the  two 
stations,  which  indicates  local  (i.e.,  longitudinal)  control  sometimes 
dominate.  The  most  likely  candidate  for  this  control  is  neutral  atmo¬ 
spheric  motion. 
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DD 

Field  Command 

Defense  Nuclear  Agcy 
ATTN:  FCP 

Field  Command 

Defense  Nuclear  Agcy 

Livermore  Br 

ATTN:  FCPRL 

Interservice  Nuclear  Wpns  Sch 
ATTN:  TTV 

Joint  Chiefs  of  Staff 

ATTN:  C3S,  Eval  Ofc 
ATTN:  C3S 

Joint  Strat  Tgt  Planning  Staff 
ATTN:  JLA 
ATTN:  JLTW-2 

National  Security  Agcy 

ATTN:  R-52,  J.  Skillman 
ATTN:  B-3,  F.  Leonard 
ATTN:  W-32,  0.  Bartlett 


DEPARTMLNT  OF  DEFENSE  (Continued) 

Linder  Secretary  of  Defense  for  Rsch  8  Engrg 
ATTN:  Strategic  &  Space  Sys  (OS) 

WWMCCS  System  Engrg  Org 
ATTN:  R.  Crawford 

DEPARTMENT  OF  THE  ARMY 

Assistant  Chief  of  Staff  for  Automation  8  Conin 

Department  of  the  Army 

ATTN:  DAAC-ZT ,  P.  Kenny 

Atmospheric  Sciences  Lab 

US  Army  Electronics  R8D  Cmd 

ATTN:  DELAS-EO,  F.  Niles 

BMD  Advanced  Technology  Ctr 

Department  of  the  Army 

ATTN:  ATC-T,  M.  Capps 
ATTN:  ATC-0,  W.  Davies 

BMD  Systems  Cmd 

Department  of  the  Army 
2  cy  ATTN:  BMDSC-HW 

Deputy  Chief  of  Staff  for  Ops  8  Plans 

Department  of  the  Army 
ATTN:  DAMO-RQC 

Harry  Diamond  Labs 

Department  of  the  Army 

ATTN:  DELHD-I-TL,  M.  Weiner 
ATTN:  Chief,  Div  20000 
ATTN:  DELHD-N-RB ,  R.  Williams 

US  Army  Chemical  Sch 

ATTN:  ATZN-CM-CS 

US  Army  Comm-Elec  Engrg  Instal  Agcy 
ATTN:  CCC-EMEO-PED,  G.  Lane 
ATTN:  CCC-CED-CCO,  W.  Neuendorf 

US  Army  Communications  Cmd 
ATTN:  CC-OPS-W 
ATTN:  CC-OPS-WR,  H.  Wilson 

US  Army  Communications  R8D  Cmd 

ATTN:  DRDCO-COM-RY ,  W.  Kesselman 

US  Army  Foreign  Science  8  Tech  Ctr 
ATTN:  DRXST-SD 

US  Army  Materiel  Dev  8  Readiness  Cmd 
ATTN:  DRCLDC,  0.  Bender 

US  Army  Missile  Intel  Agcy 
ATTN:  YSE ,  J.  Gamble 

US  Army  Nuclear  8  Chemical  Agcy 
ATTN:  Lib 
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DEPARTMENJ  OF  THE  ARMY  (Continued) 

US  Army  Satellite  Conn  Agcy 
ATTN:  Doc  Con 

US  Army  TRADOC  Sys  Analysis  Actvy 
ATTN:  ATAA-PL 
ATTN:  ATAA-TDC 
ATTN:  ATAA-TCC,  F.  Payan,  Jr 

DEPARTMENT  OF  THE  NAVY 

COMSPTEVFOR 

Department  of  the  Navy 

ATTN:  Code  605,  R.  Berg 

Joint  Cruise  Missiles  Proj  Ofc 

Department  of  the  Navy 
ATTN:  JCMG-707 

Naval  Air  Development  Ctr 

ATTN:  Code  6091,  M.  $etz 

Naval  Air  Systems  Cmd 
ATTN:  PMA  271 

Naval  Electronic  Systems  Cmd 

ATTN:  PME  117-211,  B.  Kruger 
ATTN:  PME  106-4,  S.  Kearnev 
ATTN:  PME  106-13,  T.  Griffin 
ATTN:  PME  117-2013,  G.  Burnhart 
ATTN:  Code  501A 
ATTN:  PME  117-20 
ATTN:  Code  3101,  T.  Hughes 

Naval  Intelligence  Support  Ctr 
ATTN:  N I  SC- 50 

Naval  Ocean  Systems  Ctr 

ATTN:  Code  532,  R.  Pappert 
ATTN:  Code  532,  J.  Bickel 
ATTN:  Code  5322,  M.  Paulson 
3  cy  ATTN:  Code  5323,  J.  Ferguson 

Naval  Research  Lab 

ATTN:  Code  4780,  S.  Ossakow 
ATTN:  Code  7950,  J.  Goodman 
ATTN:  Code  7550,  J.  Davis 
ATTN:  Code  4187 
ATTN:  Code  4700,  T.  Coffey 
ATTN:  Code  7500,  B.  Wald 

Naval  Space  Surveillance  Sys 
ATTN:  J.  Burton 

Naval  Surface  Wpns  Ctr 
ATTN:  Code  131 

Naval  Telecom  Cmd 

AT  1 N :  Code  341 

Office  of  Naval  Rsch 
ATTN:  Code  420 
ATTN:  Code  465 
ATTN:  Code  421 

Office  of  the  Chief  of  Naval  Ops 
ATTN:  OP  98 IN 
A 1  T N :  OP  94 ID 
AT  IN:  OP  65 


DEPARTMENT  OF  THE  NAVY  (Continued) 

Strategic  Systems  Proj  Ofc 
Department  of  the  Navy 
ATTN:  NSP-2141 
ATTN:  NSP-43 

ATTN:  NSP-2722  ,  F.  Wimberly 

DEPARTMENT  OP  THE  AIR  FORCE 

Aerospace  Defense  Cmd 
Department  of  the  Air  Force 
ATTN:  DC,  T.  Long 

Air  Force  Geophysics  Lab 

ATTN:  OPR,  H.  Gardiner 
ATTN:  OPR- 1 
ATTN:  LKB ,  K.  Champion 
ATTN:  OPR,  A.  Stair 
ATTN:  S.  Basu 
ATTN:  PHP 

ATTN:  PHI,  J.  Buchau 
ATTN:  R.  Thompson 

Air  Force  Weapons  Lab 
Air  Force  Systems  Command 
ATTN:  SUL 
ATTN:  NTYC 
ATTN:  NTN 

Air  Force  Wright  Aeronautical  Lab 
ATTN:  W.  Hunt 
ATTN:  A.  Johnson 

Air  Logistics  Command 
Department  of  the  Air  Force 
ATTN:  00-ALC/MM 

Air  University  Lib 
Department  of  the  Air  Force 
ATTN:  AUL-LSE 

Air  Weather  Service,  MAC 
Department  of  the  Air  Force 
ATTN:  DNXP,  R.  Babcock 

Assistant  Chief  of  Staff 
Studies  &  Analyses 
Department  of  the  Air  Force 

ATTN:  AF/SASC,  C.  Riqhtmeyer 
ATTN:  AF/SASC,  W.  Keaus 

Ballistic  Missile  Ofc 
Air  Force  Systems  Command 

ATTN:  ENSN,  J.  Allen 

Deputy  Chief  of  Staff 
Ops  Plans  and  Readiness 
Department  of  the  Air  Force 


ATTN: 

ATXOKS 

ATTN : 

AFXOXFD 

ATTN: 

AFXOKT 

ATTN: 

ATXOtXD 

Deputy  Chief 

of  Staff 

Rsch,  Dev,  & 

Acq 

Department  of  the  Air 

ATTN: 

AFTTUS 

ATTN: 

AFRDSP 

ATTN: 

AFRDSS 
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DEPARTMENT 

OF  THE  AIR  FORCE  (Continued) 

DEPARTMENT 

OF  ENERGY  CONTRACTORS 

Electronic 

Systems  Div 

EG&G,  Inc 

Air  force  Systems  Command 

ATTN 

J.  Colvin 

ATTN 

OCKC,  J.  Clark 

ATTN 

D.  Wright 

Electronic 

Systems  Div 

Lawrence  Livermore  Natl  Lab 

Air  Force  Systems  Conmand 

ATTN 

L-389 ,  R.  Ott 

ATTN' 

OCT-4,  J.  Deas 

ATTN 

L-31,  R.  Hager 

ATTN 

Tech  Info  Dept  Lib 

Electronic 

Systems  Div 

Air  force  Systems  Command 

Los  Alamos 

Natl  Lab 

ATTN 

YSM,  J.  Kobelski 

ATTN 

D.  Simons 

ATTN 

VSEA 

ATTN 

E.  Jones 

ATTN 

D.  Westervelt 

Foreign  Technology  Div 

ATTN 

P.  Keaton 

Air  Force  Systems  Command 

ATTN 

MS  670,  J.  Hopkins 

ATTN 

TQTD,  B.  Ballard 

ATTN 

R.  Taschek 

ATTN 

N11S,  Lib 

ATTN 

MX  664,  J.  Zinn 

Headquarters  Space  Div 

Sandia  National  Labs 

Air  Force  Systems  Command 

Livermore  Lab 

ATTN 

SKA,  D.  Bolin 

ATTN 

B.  Murphey 

ATTN 

SKY,  C.  Kennedy 

ATTN 

T.  Cook 

Headquarters  Space  Div 

Sandia  Nat 

Lab 

Air  Force  Systems  Command 

ATTN 

Org  4241 ,  T .  Wright 

ATTN 

YZJ,  W.  Mercer 

ATTN 

D.  Thornbrough 

ATTN 

Org  1250,  W.  Brown 

Headquarters  Space  Div 

ATTN 

3141 

Air  Force  Systems  Command 

ATTN 

D.  Dahlgren 

ATTN 

E.  Butt 

ATTN 

Space  Proj  Div 

Rome  Air  Development  Ctr 

DEPARTMENT 

OF  DEFENSE  CONTRACTORS 

Air  Force  Systems  Command 

ATTN 

OCS,  V.  Coyne 

Aerospace  Corp 

ATTN 

TSLD 

ATTN 

R.  Slaughter 

ATTN 

J.  Straus 

Rome  Air  Development  Ctr 

ATTN 

V.  Josephson 

Air  Force  Systems  Command 

ATTN 

I.  Garfunkel 

ATTN 

EEP 

ATTN 

D.  Olsen 

ATTN 

T.  Salmi 

Strategic  Air  Cmd 

ATTN 

N.  Stockwell 

Department  of  the  Air  Force 

ATTN 

S.  Bower 

ATTN 

DCXT 

ATTN 

DCXR,  T.  Jorgensen 

University 

of  Alaska 

ATTN 

NRT 

ATTN 

Tech  Lib 

ATTN 

XPFS 

ATTN 

N.  Brown 

ATTN 

DCX 

ATTN 

T.  Davis 

OTHER  GOVERNMENT  AGENCIES 

Analytical 

Systems  Engrg  Corp 

ATTN 

Radio  Sciences 

Central  Intel liqence  Aqcy 

ATTN 

OSWR/NED 

Analytical 

Systems  Enqrq  Corp 

ATTN 

Scty 

Department 

of  Commerce 

National  Bureau  of  Standards 

Barry  Research  Corp 

ATTN 

Sec  Ofc  for  R.  Moore 

ATTN 

J.  McLaughlin 

Department 

of  Commerce 

BDM  Corp 

National  Oceanic  &  Atmospheric  Admin 

ATTN 

L.  Jacobs 

ATTN 

R.  Grubb 

ATTN 

T.  Neighbors 

Institute  for  Telecom  Sciences 

Berkeley  Research  Assoc,  Inc 

National  Telecom  &  Info  Admin 

ATTN 

J.  Workman 

ATTN 

A.  Jean 

ATTN 

L.  Berry 

Betac 

ATTN 

W.  Utlaut 

ATTN 

J.  Hirsch 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Boeing  Co 

ATTN:  M/S  42-33,  J.  Kennedy 
ATTN:  G.  Hall 
ATTN:  S.  Tashird 

Booz-Allen  8,  Hamilton,  Inc 
ATTN:  B.  Wilkinson 

University  of  California  at  San  Diego 
ATTN:  H.  Booker 

Charles  Stark  Draper  Lab,  Inc 
ATTN:  J.  Gilmore 
ATTN:  D.  Cox 

Communications  Satellite  Corp 
ATTN:  D.  Fang 

Computer  Sciences  Corp 

ATTN:  F,  Eisenbarth 

Comsat  Labs 

ATTN:  R.  Taur 
ATTN:  G.  Hyde 

Cornell  University 

ATTN:  M.  Kelly 
ATTN:  D.  Farley,  Jr 

E-Systems,  Inc 

ATTN:  R.  Berezdivin 

Electrospace  Systems,  Inc 
ATTN:  H.  Logston 

ESL,  Inc 


ATTN: 

J.  Marshall 

General  Electric  Co 

ATTN: 

M.  Bortner 

ATTN: 

A.  Harcar 

General  Electric  Co 

ATTN: 

C.  Zierdt 

ATTN: 

A.  Steinmayer 

General  Electric  Co 

ATTN: 

F.  Reibert 

General  Electric  Co 

ATTN: 

G.  Millman 

General  Res': 

■arch  Corp 

ATTN: 

J,  Ise,  Jr 

ATTN: 

J.  Garbarir.o 

Harris  Corp 

ATTN: 

E.  Knick 

Horizons  Technology,  Inc 

ATTN: 

R,  Kruger 

HSS,  Inc 

ATTN: 

D.  Hansen 

IBM  Corp 

ATTN: 

F.  Ricci 

DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


Uni versi ty 

of  Illinois 

ATTN 

K.  Yeh 

Insti tute 

for  Defense  Analyses 

ATTN 

E.  Bauer 

ATTN 

H.  Wolfhard 

ATTN 

J.  Aein 

ATTN 

J.  Bengston 

International  Tel  &  Tel  Corp 

ATTN 

Tech  Lib 

ATTN 

G.  Wetmore 

JAYCOR 

ATTN 

J.  Sperling 

JAYCOR 

ATTN 

J.  Doncarlos 

Johns  Hopkins  University 

ATTN 

T.  Potemra 

ATTN 

J.  Phillips 

ATTN 

T.  Evans 

ATTN 

J.  Newland 

ATTN 

P.  Komiske 

Kaman  Tempo 

ATTN 

DAS  I  AC 

ATTN 

T.  Stephens 

ATTN 

W.  McNamara 

ATTN 

W.  Knapp 

Linkabit  Corp 

ATTN 

I.  Jacobs 

Litton  Systems,  Inc 

ATTN 

R.  Grasty 

Lockheed  Missiles  &  Space  Co,  Inc 

ATTN 

W.  Imhof 

ATTN 

M.  Walt 

ATTN 

R.  Johnson 

Lockheed  Missiles  &  Space  Co,  Inc 

ATTN 

Dept  60-12 

ATTN 

D.  Churchill 

ATTN 

C.  Old 

M.I.T.  Lincoln  Lab 

ATTN 

D.  Towle 

Martin  Marietta  Corp 

ATTN 

R.  Heffner 

McDonnell  Douglas  Corp 

ATTN 

N.  Harris 

ATTN 

J.  Moule 

ATTN 

W.  Olson 

ATTN 

G.  Mroz 

ATTN 

R.  Halprin 

Meteor  Comnunications  Consultants 
ATTN:  R.  Leader 
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DEPARTMENT 

OF  DEFENSE  CONTRACTORS  (Continued) 

DEPARTMENT 

OF  DEFENSE  CONTRACTORS 

Mission  Research  Corp 

Rockwell  International  Corp 

ATTN 

R.  Ki 1 b 

ATTN 

R.  Buckner 

ATTN 

Tech  Lib 

ATTN 

l(.  Hendrick 

Rockwell  International  Corp 

ATTN 

F.  Fajen 

ATTN 

S.  Quilici 

ATTN 

R.  Bogusch 

ATTN 

S.  Gutsche 

Santa  Fe  Corp 

ATTN 

D.  Sappenfield 

ATTN 

D.  Paolucci 

Mitre  Corp 

Science  Applications,  Inc 

ATTN 

G.  Harding 

ATTN 

E.  Straker 

ATTN 

C.  Callahan 

ATTN 

L.  Linson 

ATTN 

A.  Kymmel 

ATTN 

C.  Smith 

ATTN 

B.  Adams 

ATTN 

D.  Hamlin 

Mitre  Corp 

Science  Applications,  Inc 

ATTN 

M.  Horrocks 

ATTN 

S  Z 

ATTN 

W.  Foster 

ATTN 

J.  Wheeler 

Science  Appl ications,  Inc 

ATTN 

W.  Hall 

ATTN 

J.  Cockayne 

Pacific-Sierra  Rsch  Corp 

SRI  International 

ATTN 

F.  Thomas 

ATTN 

R.  Leadabrand 

ATTN 

E.  Field,  Jr 

ATTN 

D.  Neil  son 

ATTN 

H.  Brode 

ATTN 

J.  Petrickes 

ATTN 

G.  Price 

Pennsylvania  State  University 

ATTN 

R.  Tsunoda 

ATTN 

Ionospheric  Rsch  Lab 

ATTN 

A.  Burns 

ATTN 

M.  Baron 

Photometries,  Inc 

ATTN 

G.  Smith 

ATTN 

I.  Kofsky 

ATTN 

w.  Chesnut 

ATTN 

w.  Oaye 

Physical  Dynamics,  Inc 

4  cy  ATTN 

R.  Livingston 

ATTN 

E.  Fremouw 

4  cy  ATTN 

C.  Rino 

4  cy  ATTN 

J.  Owen 

Physical  Research,  Inc 

ATTN:  R.  Deliberis 

Sylvania  Systems  Group 

ATTN 

I.  Kohl  berg 

RAO  Associates 

ATTN 

R.  Steinhoff 

ATTN 

R.  Lelevier 

ATTN 

R.  Turco 

Technology 

International  Corp 

ATTN 

C.  Greifinger 

ATTN 

W.  Boquist 

ATTN 

B.  Gabbard 

ATTN 

M.  Gantsweg 

Tri-Com,  Inc 

ATTN 

W.  Karzas 

ATTN 

D.  Murray 

ATTN 

H.  Ory 

ATTN 

W.  Wright 

TRW  Devense  S  Space  Sys  Gp 

ATTN 

F.  Gilmore 

ATTN 

R.  Plebuch 

ATTN 

P,  Haas 

ATTN 

D.  Dee 

R  &  D  Associates 

Utah  State 

University 

ATTN:  B.  Yoon 

ATTN 

K.  Baker 

ATTN 

L.  Jensen 

Rand  Corp 

ATTN 

J.  Dupnik 

ATTN:  E.  Bedrozian 

ATTN 

C.  Crain 

Vi  si  dyne,  Inc 

ATTN 

C.  Humphrey 

Riverside  Research  Inst 

ATTN 

J.  Carpenter 

ATTN 

V.  Trapani 
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